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Abstract 



The miniaturization of electronic components follows a dizzying pace, and 
very soon electronic circuits will have to operate in a new regime - the meso- 
scopic regime - in which a full quantum mechanical treatment is needed. 
Furthermore, there are experimental evidences that the photosynthesis pro- 
cess fully exploit the quantum coherent transport and it is precisely this 
fact that enables natural light-harvesting systems to be so efficient. In this 
context, the system considered may have important applications in the in- 
formation technology and also in the modeling of complex systems, such as 
light-harvesting systems. 

A star graph system which consists of different chains of sites coupled with 
the origin, has been considered in the tight-binding approximation, within a 
closed and open formulation. As far as the closed system is concerned, we 
showed analytically that there are two localized states outside the normal 
Bloch band, providing the localization lengths as a function of the system 
parameters. Regarding the open system, we used an approach based on 
the effective non-Hermitian Hamiltonian in order to describe the coupling 
between the internal states and the environment. In this framework, we 
studied the eigenvalues and the localization lengths of the eigenstates by 
changing the coupling with the external world. We found that the degree 
of opening weakly affects the localized states and we also gave an analitical 
expression for the decay rates in the limit of small coupling. 

On the other hand, transport properties are very sensitive to the degree of 
opening of the system: an analytical estimation of the value of the coupling 
at which the superradiance transition occurs, which is also valid for chains 
of different lengths among them and with different couplings with the origin 
has been found. Furthermore, we have shown by numerical simulations that 
(for large number of sites in each chain) the maximum of the integrated 
transmission coincides with the value of the coupling at which the transition 
to superradiance occurs. 

This thesis is organised as follows: in Chap.l we provide an introduction 
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about the nanoworld and the importance of the system that we have studied; 
in Chap. 2 we present a derivation of the effective Hamiltonian. In Chap. 3 
we introduce the star graph model and the solution of this system for the 
closed case (no external coupling). In Chap. 4 we study in detail the open 
system, with particular attention to transport properties and superradiance 
transition. In Chap. 5 we present a preliminary result about the disorder and 
in Chap. 6 we summarize the main results obtained in this thesis. 
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'One day sir, you may tax it. " 




Michael Faraday's reply to British 
politician, William Gladstone, when asked 
of the practical value of electricity (1850) 



Introduction 



1.1 The Nanoworld revolution 



Nanotechnology is the study of manipulating matter on an atomic and molec- 
ular scale. Generally, nanotechnology deals with structures sized between 1 
to 100 nanometre in at least one dimension, and involves developing ma- 
terials or devices possessing at least one dimension within that size. Thus, 
nanometer-scale materials straddle the border between the molecular and the 
macroscopic. They are small enough to exhibit quantum properties reminis- 
cent of molecules, but large enough for their size and shape to be designed 
and controlled. The potential of nanoscale materials is almost limitless, but 
scientists must first overcome two fundamental challenges. The first is phys- 
ical: how can one control individual atoms in a nanosolid and then assemble 
them into real- world system? The second is conceptual: how does one at- 
tack problems too big to be solved by brute force calculation but too small to 
be tackled by statistical method? Overcome these problems will enable for 
nanoscience and nanotechnology to revolutionize science and technology in 
ways that will make the world 50 years from now unrecognizable compared 
with today. In addition, the nanotechnology revolution meets three very spe- 
cific scientific needs, from three different branches of science. The first is the 
rapid advances in molecular biology that have completely changed people's 
understanding of life over the past 30 years; the second is the evolution of 
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chemistry from the study of single atoms and molecules to the fabrication 
of very large complexes such as quantum dots and proteins. The third, and 
perhaps the most important need comes from Information Technology (IT): 
the continuous request to be able to handle increasing volumes of informa- 
tion imposed a rapid miniaturization of electronic components, which began 
with the birth of the transistor and the microchip. Already in 1965, Gor- 
don Moore, co-founder and president of Intel, while preparing a speech for a 
meeting, noticed that the number of transistors that can be placed inexpen- 
sively on an integrated circuit doubles approximately every two years: this is 
the well-know Moore's law, that describes a long-term trend in the history of 
computing hardware. This trend has continued for more than half a century 
and is expected to continue until 2015 or 2020 or later. This process is really 
amazing: there are no other branches of industry where advances in tech- 
nology follow an exponential process, and for so long time. Extrapolation 
of Moore's law suggests that in the next 20-30 years, electronic circuit ele- 
ments will shrink to the size of single atoms. Even before this fundamental 
limit is reached, electronic circuit will have to operate in a new regime in 
which quantum mechanics cannot be ignored. New devices that store, pro- 
cess and communicate information in a faster way will have to be invented 
in order to extend the IT revolution. No matter what device concepts are 
pursued, smaller devices must eventually embrace a profound change in the 
way we think about computing. Present devices average the behaviour of a 
great number of quantum particles. Make any of these devices small enough 
and there will be only a few quantum particle in the system. In this limit, 
the laws of Quantum Mechanics manifest themselves most vividly. Thus, 
researchers need to traverse the threshold and look beyond the hmits of clas- 
sical physics. And it is in this fervent and extremely important context that 
the physics of mesoscopic systems - systems in the intermediate size range 
between microscopic and macroscopic - moves its footsteps. 



7 



1.2. MESOSCOPIC PHYSICS 




Figure 1.1: The Moore's law and the development of Intel processors. The shrinking of 
devices is impressive; the electronic devices have nearly reached molecular scale (molecules 
commonly used as building blocks for organic synthesis have a dimension of a few A to 
several dozen A) and will achieve very quickly lengths comparable to the atomic scale 
(Bohr radius of hydrogen has approximately a value of 0.53 A). At present, everyone can 
buy a PC with a processor built with 32nm (320 A) manufacturing process technology. 
Figure taken from the IDF - Intel Developer Forum - which took place from September 
22 through September 24th, 2009. 

1.2 Mesoscopic physics 

Mesoscopic physics is a rather young branch of science. It started about 25 
years ago and it enjoys the unique combination of being able to deal with 
and provide answers to fundamental questions of physics while being relevant 
for applications in the not-too distant future, as explained in the prevoius 
section. In fact, some of the experimental possibilities in this field have been 
developed with an eye to reducing the size of electronic components. 
More formally, mesoscopic physics is a sub-discipline of condensed matter 
physics which deals with materials of an intermediate length scale. The 
length scale of such materials can be put between the size of few atoms (or 
a molecule) and of materials measuring microns. The lower limit can be also 
defined as being the size of individual atoms. At the micron level are bulk ma- 
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terials. Mesoscopic and macroscopic objects have in common that they both 
contain a large number of atoms. Whereas average properties derived from 
its constituent materials describe macroscopic objects, as they usually obey 
the laws of Classical Mechanics, a mesoscopic object, by contrast, is affected 
by fluctuations around the average, and is subject to Quantum Mechanics. 
Thus, we are in the so-called quantum realm, where quantum mechanical ef- 
fects become important and not negligible. Typically, this means distances 
of 100 nanometers or less. Even if this distance seems extremely small, it 
can be achieved quite easily today. Indeed, modern lithographic techniques 
allow researchers to pattern materials into devices with dimensions down to 
approximately 30 nanometers, roughly 100 atoms across. These systems are 
so small that they can be easily simulated with a personal computer, as we 
will do in Chapter [3] and Chapter |4j 




Figure 1.2: Electron microscope pictures of nanoparticle chains. Si nanoparticle chain 
(left) and Fe nanoparticle chain (right). Figures taken from Ref.[5j and Ref. respec- 
tively. 

1.2.1 Mesoscopic conductors 

While the resistance of an electrical element measures its opposition to the 
passage of an electric current, the electrical conductance measures how easily 
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electricity flows along a certain path (it is the inverse of the electrical resis- 
tance). Applying a voltage V across the point contact of a circuit induces a 
current to flow, the magnitude of this current is given by / = GV, where G 
is the conductance of the contact. It is well-know that the conductance of a 
conductor is directly proportional to its cross sectional area W and inversely 
proportional to its length L; namely, 

where the conductivity a is a material property of the sample independent 
of its dimensions. How small can we make, W and/or L, before the ohmic 
behavior breaks down? This question led to important developments, both 
theoretical and experimental, in our understanding of the meaning of re- 
sistance at the microscopic level. Small conductors whose dimensions are 
indermediate between the microscopic and the macroscopic are called meso- 
scopic. They are much larger than microscopic objects like atoms, but not 
large enough to be "ohmic'Q For this type of conductors Equation ( 1.1 ) is no 
longer valid and two important points should be considered. Firstly, there 
is an interface resistance independent of the length L of the sample. Sec- 
ondly, the conductance does not increase linearly with the width W. Instead 
it depends on the number of transverse modes M in the conductor and goes 
down in discrete steps. These corrections lead to the famous Laundauer's 
formul(r\ 

G = -^MT (1.2) 

The factor T represents the average probability that an electron injected at 
one end of the conductor will transmit to the other end. Now, a question 
arises: when a conductor shows an ohmic behavior and when not? Do these 
characteristic length scales, enable us to roughly define a borderline between 
the two regimes? 



-'^A notable example of this type of conductors is a Quantum Point Contact (QPC), a 
narrow constriction between two wide electrically-conducting regions, of a width compa- 
rable to the electronic wavelength (from nanometer to micrometer), see Fig. (1.3). 

formal derivation of the Laundauer's formula can be found in Refs. [Sj and [1]. 
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Figure 1.3: This image shows the coherent flow of electrons through a Quantum Point 
Contact (QPC) formed in a two dimensional electron gas (2DEG) inside a GaAs/AlGaAs 
heterostructure. The gray areas outline the position of electrostatic gates. Scanned probe 
microscope images on the outside agree well with theoretical simulations inside. The 
fringes are spaced by half the electron wavelength. Source: Robert Westervelt and Eric 
Heller at Harvard University, Arthur C. Gossard at UC Santa Barbara, Physics Today 
56(12) (2003). 

1.2.2 Quantum coherent transport 

A conductor usually shows ohmic behavior if its dimensions are much larger 
than each of three characteristic length scales: 

1. the de Broglie wavelength, which is related to the kinetic energy of the 
electrons; 

2. the inelastic mean free path, which is the distance that an electron 
travels between two collisions; 

3. the phase-relaxation length, which is the distance that an electron trav- 
els before its initial phase is destroyed. 

These length scales change widely from one material to another and are 
also strongly affected by temperature, magnetic fields etc. For this reason, 
mesoscopic transport phenomena have been observed in conductors having 
a wide range of dimensions from a few nanometers to hundreds of microns. 
Now, we want to provide some realistic estimations for the characteristic 
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lengths introduced above. In our examples, we will consider a 2D system 
because this is a situation that often occurs in experiments; for instance, 
quantum point contacts are formed in 2-dimensional electron gases (2DEG), 
e.g. in GaAs/AlGaAs heterostructure^ 

De Broglie wavelength (A) 

It is possible to sho"«|^that the Fermi wavenumber kf goes up as the square 
root of the electon density. The corresponding wavelength goes down as the 
square root of the electron density n^: 

Xf = — = J— 1.3 
kf \ n. 

For an electron density of 5 ■ 10^^ /cm^, the Fermi wavelength is about 35 
nm, comparable with characteristic lengths achieved in the industrial pro- 



duction of microprocessors, see Sec. (1.1). At low temperatures the current 



is carried mainly by electrons having an energy close to the Fermi energy so 
that the Fermi wavelength is the relevant length. Other electron with less 
kinetic energy have longer wavelengths but they do not contribute to the 
conductance. 



Inelastic mean free path (Lm) 

An electron in a perfect crystal moves as if it were in vacuum but with a dif- 
ferent mass (called effective mass). Any deviation from perfect crystallinity 
such as impurities, lattice vibrations (phonons) or other electrons leads to 
"collisions" that scatter the electron from one state to another thereby chang- 
ing its momentum. The mean free path, Lm, is the distance that an electron 
travels before its initial momentum is destroyed; that is, 

Lm = VfTm (1.4) 

■^See for example Ref.[6j. 
^See Ref.[3] for more details. 
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where is the momentum relaxation time and Vf is the Fermi velocity. The 
Fermi velocity is given by 

^ % ^ A^2^ = 3 • lO^cm/s if = 5 • lO^Vcm^ (1-5) 

mm 

Assuming a momentum relaxation time lOOps we obtain a mean free path 
Phase-rel£Lxation length(L<^) 

The most relevant length scale quantifying a mesoscopic system is the phase 
coherence length L^, the length scale over which the carriers preserve their 
phase information. This phase coherence length is, on the other hand, highly 
sensitive to temperature and sharply decreases with the rise of temperature. 
One way to visualize the distruction of phase is in term of an ideal experiment 
involving interference. For example, suppose we split a beam of electrons into 
two paths and then recombine them. In a perfect crystal, the two paths would 
be identical resulting in constructive interference. However, if we introduce 
impurities and defects randomly into each arm, the two paths arc then no 
longer identical and the interference may be destructive. Moreover, there arc 
other important situations that affect the phase-relaxation time such as 
the effect of a dynamic scatter like lattice vibrations (phonons) or electron- 
electron interactions. Therefore, to be in the mesoscopic regime we need 
low temperatur (on the order of liquid He, ~ AK) in order to neglect phase 
randomization processes caused by phonons. 

In this thesis, we consider the coherent quantum transport, which means 
that the results obtained have to be apphed to samples of a length L such 
that the electron's phase coherence length (the typical distance the elec- 
tron travels without losing phase coherence) is larger than or comparable to 
the system's size L. Phase coherence is affected by the couphng of the elec- 
tron to its environment, and phase-breaking processes involve a change in 
the state of the environment. In most cases phase coherence is lost in inelas- 
tic scatterings, e.g., with other electrons or phonons, but spin-flip scattering 
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from magnetic impurities can also contribute to phase decoherence. Elastic 
scatterings of the electron, e.g., from impurities, usually preserve phase co- 
herence and are characterized by the elastic mean free path. increases 
rapidly with decreasing temperature, and for L ~ Ifim, an open system typ- 
ically becomes mesoscopic below ~ lOOmK. In a mesoscopic sample, the 
description of transport in terms of local conductivity breaks down, and the 
whole sample must be treated as a single, coherent object]^ 
Thus, systems satisfying the condition L have to be treated quan- 

tum mechanically, in contrast with the macroscopic objects where usually 
the laws of classical mechanics are used. We want to point out that these 
systems are not so far from the size already achieved by the manufactur- 
ers of electronic components; since 2009 everyone can buy a personal com- 
puter equipped with microprocessors with 32-nm technologj]^ Therefore, 
the assumption L <^ L^, which are at the basis of our study of coherent 
quantum transport is fully justified and are not as far as can be imagined 
at first glance by important applications in everyday life. Furthermore, sev- 
eral spectacular effects appear as a consequence of quantum phase coherence 
of the electronic wave functions in mesoscopic systems like one-dimensional 
(ID) quantum wires, quantum dots where electrons are fully confined, two- 
dimensional (2D) electron gases in heterostructures, etc. For example, we 
may mention Aharonov-Bohm effect. Universal Conductance Fluctuations, 
Persistent Current, Anderson Localization and a new and very interesting 
phenomenon which is the coherent quantum transport in photosynthesis. Let 
us analyze this last phenomenon in more detaiQ 

^See Ref.[i^ from more details. 

^Moreover, in September 2009, Intel disclosed that it had the world's first working 
22nm silicon technology. Specifically, Intel disclosed that it had fimctional SRAM test 
chips with 364 Mbits, a whopping 2.9 billion transistors per chip. These SRAMs had the 
smallest cell in working circuits to date, at 0.092 square microns. The 22nm process is on 
track for production in the 2"^^ half of 2011, two years after start of 32nm high volume 
production. 

''See Ref.[S] and the references therein. 



14 



1.2. MESOSCOPIC PHYSICS 




Figure 1.4: Like ripples on water, energy absorbed from sunlight moves through the 
photosynthetic complex (gray-and-green molecular structure) with a wavelike motion. In 
this artist's rendering, that motion creates interference signals that are analyzed in a 
spectrum (left). (Image of Greg Engel, Lawrence Berkeley National Laboratory, Physical 
Biociences Division) 



Coherent quantum transport and photosynthesis 

In organisms ranging from blue algae to giant sequoias, complicated assem- 
blies of molecules of the pigment chlorophyll absorb sunlight's photons and 
channel their energy to enable the plants to turn water and carbon dioxide 
into oxygen and sugars. The efficiency of photosynthesis, as this process is 
called, has long astounded scientists. Virtually every photon absorbed by 
chlorophyll initiates a photosynthetic reaction. Plants use up to 95 percent 
of the light that strikes them, whereas commercial solar panels use less than 
30 percent. But how do plants achieve this amazing result? 
Before answering this crucial question, let us look at the mechanism of pho- 
tosynthesis. Photosynthesis is initiated by the excitation, through incident 
light, of electron in pigment molecules - chromophores - such as chlorophyll. 
This electronic excitation moves downhill from energy level to energy level 
through the chromophores before being trapped in a reaction center, where 
its remaining energy is used to initiate the production of energy-rich car- 
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bohydrates. The mechanism of energy transfer! through chromophore com- 
plexes has generally been assumed to involve incoherent hopping, thus doing 
a random walk with a general downhill direction. But at each hopping, the 
excitation might dissipate as waste heat, so scientists did not understand how 
the process could be so efficient. The solution was experimentally discovered 
by Engel et al^ they found that groups of chlorophyll molecules spend a 
surprisingly long time in a superposition of states. In the experiment, the 
team froze (77 kelvin) chlorophyll complexes from blue algae and shot them 
with sequences of ultrashort laser pulses, each lasting just 40 femtoseconds. 
Three pulses excited the molecules, and a fourth pulse detected interference 
patterns. The complexes stayed in a superposition of states for more than 
600 femtoseconds after receiving the pulses. In other words, the electronic 
excitation that transfers the energy downhill does not simply hop incoher- 
ently from state to state, but samples two or more states simultaneously. 
Thus, the photosynthesis process fully exploit (also at room temperature!) 
the quantum coherent transport and it is precisely this fact that enables nat- 
ural light-harvesting systems to be so efficient. Concluding, we can easily 
understand how the quantum coherent transport could provide a clean so- 
lution to mankind's energy requirements. Nature already knows to do that, 
now it's up to us. 

1.3 The importance of the star graph system 

The miniaturization of electronic components follows an impressive pace, and 
very soon electronic circuits will have to operate in a new regime - the meso- 
scopic regime - in which a full quantum mechanical treatment is needed, in 
contrast to the macroscopic objects where usually the laws of classical me- 
chanics are used. Furthermore, there are experimental evidences that the 
photosynthesis process fully exploit the quantum coherent transport and it 
is precisely this fact that enables natural light-harvesting systems to be so 

*See Ref.|9] for more details. 
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Figure 1.5: The star graph system. 

efficient. In this context, the system considered may have important apph- 
cations in the information technology and also in the modeling of complex 
systems, such as light-harvesting systems. 

Up to now, in the literature were considered only very simple systems, such 
as Id, 2d and 3d chains; in this thesis, we go beyond because we consider a 
more complex system, but still analytically treatable, which is closest to the 
complex systems existing in nature than chains considered so far. 
More precisely, in this thesis, we will study in details what we named star 
graph system, a system which consist of h chains of sites with a common ver- 
tex at the point O, the origin. The only coupling between the chains comes 
through the origin. Furthermore, when we will considerer an open system. 



the last site of each chain will be coupled to the external world. In Fig. (1.5) 



we show a representation of this system. In this thesis, we will present the 
results both for the closed and the open system; we emphasize the fact that 
all theoretical and numerical results related to the above system contained 
in this thesis are brand new and unpublished. In fact, it is the first time that 
the star graph system has been studied in such detail; up to now, such a 
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detailed study of a similar system had been done only for a one-dimensional 
chain, see Ref.[I]. Now, we want to share with the reader the importance of 
the star graph system, and also the coherent quantum transport, in such a 
system. Here, we present some crucial points: 

• in the mesoscopic regime, electronic transport cannot be investigated 
by using the conventional Boltzmann transport equation since at this 
length scale quantum phase coherence plays an important role and a 
full quantum mechanical treatment is needed. 

• The rapid miniaturization of electronic devices will require understand- 



ing the mechanism of coherent quantum transport, see Sec. (1.1) and 



Sec. (1.2); furthermore, if we will build a quantum computer, we will 
have to know in detail the quantum coherent transport, in order to fully 
exploit the potential of quantum computing such as the quantum paral- 
lelism. In this context, the star graph system could be used as quantum 
wires to connect different components of a quantum computer. 

The notion of a closed physical system is always an idealization be- 
cause the interaction with the outside world, including the measuring 
apparatus, cannot be removed completely; furthermore, the transport 
properties depend strongly on the degree of openess of the system. 
Therefore, our studies on the star graph open system are of great im- 
portance in understanding the interaction of our quantum system with 
the environment. The effective non-Hermitian Hamiltonian approach 
to open quantum system used in this thesis has been shown to be a very 
effective tool to control the effect of the opening of a quantum system. 
A clear example of non trivial effect of openess on transport properties 
is the phenomenon of superradianc^ recentely shown to occurs in a 
paradigmatic model of quantum coherent transporlj^ 



^See more detail in Chapj2j Sec.(E4| 
lOSee Refs.l] and g]. 
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• In this thesis, we will demonstrate the existence of localized states in 
branching regions; we will show that these states are very little affected 
by the opening of the system and thus have a rather large decay time, 



feature of the star graph system, could be used as memory cells. 

Concluding, we can state that the star graph system studied has lots of new 
items that could be utilized in many different areas, such as condensed mat- 
ter, nanotechnology, molecular biology, information technology and quantum 
computation. 



Figure 1.6: (a) Optical micrograph showing Cr/Au electrical leads to five crossed SWNT 
(single-walled carbon nanotubes) devices, (b) Tapping mode AFM image (amplitude 
signal) of the area indicated by the white box in (a). Figures and caption taken fi'om 




Accordingly, these localized states, distinctive 




Ref.[IO]. 
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There must be no barriers for freedom of 
inquiry. There is no place for dogma in sci- 
ence. The scientist is free, and must be free 
to ask any question, to doubt any assertion, 
to seek for any evidence, to correct any er- 
rors. 



Robert Oppenheimer - American physicist 




Open quantum systems: 
the effective Hamiltonian approach 



Transport properties depend strongly on the degree of openess of the system. 
In important apphcations, the effect of the opening is large and cannot be 
treated perturbatively. In a typical situation, we have a discrete quantum 
system coupled to an external environment characterized by a continuum 
of states. Elimination of the continuum leads to an effective non-Hermitian 
Hamiltonian. Analysis of the complex eigenvalues of the effective Hamilto- 
nian reveals a general phenomenon, namely, the segregation of decay widths 
(corresponding to the imaginary parts of the complex eigenvalues). Generi- 
cally, at weak coupling, all internal states are similarly affected by the opening 
and acquire small decay widths, resulting in narrow transmission resonances. 
As the coupling increases and reaches a critical value, the resonances overlap, 
and a sharp restructuring of the system occurs. Beyond this critical value, a 
few resonances become short-lived states, leaving all other {long-lived) states 
effectively decoupled from the environment. This general phenomenon is 
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referred to as the super-radiance transition, due to its analogy with Dicke 
super-radiance in quantum optic^ In this chapter, we will present a deriva- 
tion of the general form of the effective non-Hermitian Hamiltonian starting 
from standard scattering theory. 



2.1 The Lippmann-Schwinger equation 

We assume that the Hamiltonian can be written as 

H = Ho + V (2.1) 
where Hq stands for the kinetic-operator 

In the absence of a scatterer, V would be zero, and an energy eigenstate 
would be just a free particle state The presence of V causes the energy 
eigenstate to be different from a free-particle state. However, in presence 
of elastic scattering we are interested in a solution of the fuU-Hamiltonian 
Schrodinger equation with the same energy eigenvalue. The free particle 
state 10) is the energy eigenket of Hq: 

HM) = E\ct>) (2.3) 

where E is the energy of the system. The basic Schrodinger equation to be 
solved is 

{H, + V) 1^) = (2.4) 



We look for a solution to Eq.(2.4) such that, as — i- 0, we have 



where 10) is the solution to the free-particle Schrodinger equation (2.3) with 
the same energy eigenvalue. 



The solution of Eq.(2.4) is: 



I^W) = 10) + ^ I . .^ V-I^W) (2.5) 



iSee Ref.[Il]. 
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where e is a small real parameter. This is known as the Lippmann-Schwinger 
equation, named after Bernard A. Lippmann and Julian Schwinger. 
Now we introduce the Green's functions of the system. Defining z = E + iCQ, 
the function 

= ^ (2.6) 

is called Green's function (or resolvent, or propagator) of H and fulfills the 
equation 

{E-H)G = 1 (2.7) 

Similarly, the function 

9i^) = V^ (2.8) 
is the Green's function of the operator Hq and satisfies 

{E-Ho)g = l (2.9) 

Now, we want to find out a relation between the two Green's functions. 
Consider the following equation 

g-\z) - G-\z) = {z- Ho) -{z-H) = -Ho + H = V. (2.10) 



Multiplying Eq.(2.10) to the left by g{z) and to the right by G{z) gives 

g{z) [g-\z) - G-\z)] G{z) = G{z) - g{z) = g{z)VG{z), (2.11) 
from which follows 

Giz) = giz) + giz)VG{z) (2.12) 

This is the Dyson equation, an equation giving the relation between the two 
propagators g{E) and G{E). 



^ This is a shortened form which stands for 



z — hm {E + ie) 

£^0 + 
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On the other hand, multiplying Eq.(2.10) to the left by G{z) and the right 
by g{z) we obtain 

Giz)=giz) + Giz)Vgiz) (2.13) 



which is another formulation of the Dyson equation, equivalent to Eq.(2.12). 



Coming back to Eq.(2.12), multiplying to the right by V we obtain 

GV = gV + GVgV 



which can be written as 



'l + GV){l-gV) = l 



Now, the Lippman-Schwinger equation (2.5) reads 



(l-(7\/)|V^W) 



(2.14) 



(2.15) 



(2.16) 



that, multiplying to the left by (1 + GV) and taking into account Eq.(2.15), 
gives 

|^W) = (l + Gy)|0) (2.17) 
It is convenient to introduce a matrix T such that 



(2.18) 



where \(f>a,b) and \4'a,b) stand for the eigenkets of the system in the state a 
and b, respectively. T will be called the transition matrix, and the element 
T""^ the amplitude for the transition a ^ b. 



Multiplying Eq.( pJ7| ) to the left by V we get 

y|^(+)) = (V + VGV) 



that is, using the definition of the transition matrix (2.18), 



,|r|0a) = {V + VGV) \4>a) 



that implies 



r = v + VGV = v + v 



1 

E-H + iO 



-V 



(2.19) 



(2.20) 



(2.21) 
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2.2 The effective Hamiltonian 



Consider a discrete quantum system described by intrinsic basis states cou- 
pled to a continuum of states. First of all, let us divide the Hilbert space of 
the full Hermitian Hamiltonian H of the system into two mutually orthogo- 
nal subspace Sp and Sq with the use of the projection operators P and 
respectively. The subspace Sp involves the internal states {|^)}; n = 1..N 
and the subspace Sq is composed by the continuum channel states {\c; E)} 
where c = 1..M is a discrete quantum number labeling M channels and E is 
a continuum quantum number representing the energy. 
Keeping in mind the division of the Hilbert space of H, the operator H can 
be written in the form 

Hpp HpQ 
Hqp Hqq 

where we used the notation Hxy = XHY . Now, splitting the operator if, 
we obtain 



H 



(2.22) 



H = Ho + V 




(2.23) 




Application of the Dyson equation (|2.12|) gives 

Gf _ ^ 

V gQQ J 

^( 9PP \ 

\ Qqq J \ Hqp J \ Gqp Gqq 

(2.24) 

where gpp{z) = [z - Hpp)'^ and qqq = {z - HqQy^. Performing the 
algebrical matrix operations indicated in the above expression, it follows 




Gpp{z) = gpp{z) + gpp{z)HpQGQp{z) 
Gqp{z) = gQQ{z)HQpGpp{z) 



(2.25) 
(2.26) 
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(and similar expressions for exchange of the subscripts P and Q). With 

straightforward properties of matrix multiphcation we have 

1 



Gpp{E) 



E-H 



pp 



H 



1 



(2.27) 



PQ- 



QP 



E-Hqq + 

Until now we have performed exact algebraic transformations. 



The interpretation of equation (2.27) is immediate. The Green's function in 



the space Sp is determined by an effective Hamiltonian which is given by 

(2.28) 



n{E) = Hpp + HpQ 



1 



E-Hr 



iO 



-H, 



QP 



The physical meaning of the effective Hamiltonian is self-explanatory: besides 
Hpp, it contains the effect of an excursion from Sp to Sq, a propagation 
within Sq, and an excursion back to Sp. 

Suppose that the splitting of the original space spanned by H is such that 
in the subspace Sq the Green's function is known. We can eliminate (or as 
it is commonly said in the literature decimate) all the states of the subspace 
Sq, considering within the subspace Sp the effective Hamiltonian 

n{E) = Hpp + ^pp{E) (2.29) 

with 



T.pp{E) = HpQ- 



1 



:H 



The operator 'Lpp{E), whose origin is linked to the elimination of the sub- 
space Sq, is called self-energy operator; the self-energy operator T,pp{E), 
added to Hpp, gives the effective Hamiltonian on the preserved subspace Sp, 
now formally decoupled from the subspace Sq. 



Coming back to Eq.(2.28), we note that the effective Hamiltonian is in term 



of the projection operators. These are the explicit structures of the projection 
operators 



TV 



(n\ 



P = J]|n)(r, 

n=l 

Q = J2j dE'\c,E'){c,E') 



(2.31) 
(2.32) 
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so that 

P + Q = l (2.33) 

and 

PQ = QP = (2.34) 

In order to put Eq.(2.28) in a more transparent form we substitute into 
it the definitions of the projection operators and, after a straightforward 
calculation, we obtain 

C T,S 

Multiplying from the left by (m| and from the right by \n) which are two 
internal states, we get 

{mrH{E)\n) = {m\H\n) + j ^^'^^^T^^ ^^'^^^ 

c 

where we have defined 

A'l^{E') = {n\H\c,E') (2.37) 
which is the transition amplitude between the intrinsic states and the contin- 



uum. The integral in Eq.(2.36) can be further decomposed into its Hermitian 
part (principal value) and the remaining-non Hermitian part by using the 
Sokhotski-Plemelj formula]^ 

^ J E-E' + iO ^ J E-E' 

(2.38) 

-^n J2 AUE)Al{Er 

c(open) 

where the notation VV refers to the Cauchy principal value. Then, the 
effective Hamiltonian for the intrinsic system, which fully takes into account 



Let / be a complex-valued function and let a and b be real constants with a < < 6. 
Then ^ ^ 

where VV denotes the Cauchy principal value. 
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its opening to the outside, can be written as 

H{E) = Hpp + A{E) - -W{E) (2.39) 

with 

Wnm = 2n J2 Al{E)AUEr (2.40) 

c(open) 

where the sum is hmited to the open channels, and 

A.. = E / dE'^^^^^^^ (2.41) 

Assuming A {E) and W (E) are smooth function of the energy, their en- 
ergy dependence can be neglected if the region of interest is concentrated 
in a small energy window. Therefore, the amplitudes can be taken as 
energy-independent parameters. Taking into account all these considera- 



tions, Eq.(2.39) becomes 

M 

U = Hpp- -W Wmn = (2-42) 

c=l 

Under time-reversal invariance, both Hpp and W are real symmetric matri- 
ces, then the coupling amplitudes between intrinsic states \n) and chan- 
nels c can be taken as real. 



2.3 The scattering matrix 



The effective Hamiltonian (2.39) determines the scattering matrix of the 



system. First of all, we rewrite Eq.(2.23) in a different, but equivalent way: 



H = H 



pp 



H 



QQ 



H 



PQ 



H 



QP 



(2.43) 



where we used the notation Hxy = XHY previously introduced. 
The first two terms Hy^ = Hpp + Hqq are the part of the full Hamiltonian 
H acting within the respective subspaces P and Q. Similarly, the last two 
terms V = Hpq + Hqp act across the two subspaces. Our full Hamiltonian 
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can now be written as H = H^j + V . 

From standard scattering theory, we know that the scattering matrix is de- 
fined as 

S = l-iT (2.44) 



where T is the transition matrix aheady defined in Eq.(2.18). 



Using Eq.(2.21), we obtain the following relation for the transition matrix: 

1 



T=V + VGV = V + V- 



-V 



E-H + iO 

where G is the Green's function of the full Hamiltonian H. 
The transition amplitude for the process from channel a to 6 is 

T"' (E) ={b,E\r\a,E) 

= {b, E\V\a, E) + (6, E\VGV\a, E) 



(2.45) 



(2.46) 



Now, we insert V = HpQ+Hgp and P = J2n i^^^ above expression 

and by exploiting the orthogonality of the two subspace Sp and Sq, we find 



T^'iE) = {b,E\HQpPGPHpQ\a,E) 



(2.47) 



We note that in the equation above the propagator PGP appears which is 
equal to Gpp. The latter can be seen as a projection of the full propagator 
G into the internal subspace, after the elimination of channel variables. 



Using both Eq.(2.27) and (2.28), the equation (2.47) reads 



N 

T''\E) = ^ {h,E\H\n){n\Gpp\m){m\H\a,E) 

n,m=l 

= V Ai* (E) ( ^——] Al, (E) 

n,m=l ^ ^ / / nm 



(2.48) 



The propagator {E — Ti) ^ in the scattering amplitude T"'' does not depend 



on a specific reaction and contains the full effective Hamiltonian (2.39) with 



the same amplitudes A'^ as those determining the entrance and exit channel 



in Eq.(2.48). This guarantees the unitarity of the iS-matrix since the virtual 



processes of evolution of the open system to and from the continuum channels 
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are included in all orders in the propagator. From the transition amplitude 
T, we define the transmission T"-^{E) which gives us the transmission from 
channel a to b: 

J^\E) = \r''\E)\'^ (2.49) 

Now, we want to write T"''' in a different way in order to show that the 
eigenvalues of "H coincide with the poles of the scattering matrix 5 = 1— iT. 
First of all, we diagonalize the effective non-Hermitian Hamiltonian "H. Its 
eigenfunctions |r) and (f| form a bi-orthogonal complete set 

n\r) = Sr\r) {r\n = {r\S; (2.50) 

and its eigenvalues are complex energies 

£r = Er - ^Tr (2.51) 

corresponding to resonances centered at Ej. with widths which determine 
the lifetime of the resonances, ~ h/Tr- The decay amplitudes Af^^ are 
transformed according to 

y^- = J2'^n{n\r) = J](r1m)A^ (2.52) 

n m 

and the transmission amplitudes are given by 

^ 1 

r'^HE) = J2^r^^^'r (2.53) 

■r=l ' 

The complex eigenvalues Sr of "H thus coincide with the pole of the scattering 
matrix S. From this consideration, it is clear that the properties of the 
complex eigenvalues of the effective Hamiltonian Ti are extremely important 
for understanding the transport properties of the system. 



2.4 Transition to superradiance 



In order to understand what is the transition to superradiance, we can con- 

(2.54) 



sider a simplified version of Eq.(2.42): 



'H = Ho- -W 



29 
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where 

W,, = J2 "^"^ (2-55) 

c{open) 

and where 7 is a parameter that controls the couphng strength with the ex- 
ternal world (which we assume to be of the same order of magnitude for all 
the intrinsic states), and the basis states \i) are chosen to be the eigenstates 
of Hq, with eigenvalues Eq. 

As long as 7 is small (7 <^ 1), the second term of the Hamiltonian % in 



Eq.(2.54) can be considered as a small perturbation of Hq. This condition is 
always fulfilled if the average width (F) is much smaller than the average dis- 
tance D between neighboring resonance states. In this case, the nondiagonal 
matrix elements of "H are small and the individual resonances are isolated. 
So, we have that the first-order complex eigenvalues of % are 

£, = El- IpWu (2.56) 

In the opposite case of large 7 (7 ^ 1), the matrix W determines the behavior 
of the system. Hq can be viewed as a perturbation acting on W. From 



Eq.(2.55) it is evident that the rank of W is M, the number of open channels, 
and, hence, M is also the rank of l-L. This implies that W has only M nonzero 
eigenvalues for M < N. Thus, only M states will have a decay width in the 
limit of large coupling, while all others will have zero width to first order. 
Therefore, as the coupling increases, all widths initially increase linearly with 
7, but at large coupling only M of the widths continue to increase, while 
the remaining N — M widths approach zero. The N — M states almost 
decoupled from the continuum of decay channels become long-lived (trapped) 
while M states take almost the whole coupling strength and become short- 
lived (super-radiant). This general phenomenon is referred to as the super- 
radiance transition, due to its analogy with Dicke super-radiance in quantum 
optics, see Ref . |14] . Therefore, it is clear that a transition between these two 
regimes may take place at a critical value of 7. Roughly, the transition occurs 
when 

7/L> ^ 1 (2.57) 
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where D is the mean level spacing of Hq. Note that the qualitative criterion 
j/D ^ 1 for the transition to superradiance is valid in the case of uniform 
density of states and negligible energy shift; when the density of states is not 
uniform, the transition to superradiance occurs as a hierarchical process, see 
Ref. [H]. 
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A good idea has a way of becoming simpler 
and solving problems other than for which 
it was intended. 

Robert Tarjan - computer scientist 

Star graph: 
the closed tight-binding model 




3.1 Overview of tight-binding method 

The tight binding method., suggested by Bloch in 1928, is an approach to the 
calculation of electronic band structure using an approximate set of wave 
functions based upon superposition of wave functions for isolated atoms lo- 
cated at each atomic site; these wave functions can be physically interpreted 
as atomic orbitals. Orbitals of neighboring (or not neighboring) sites are 
connected by what is referred to as a hopping matrix element or an over- 
lap integral. The tight-binding approximation deals with the case in which 
the overlap of atomic wave functions is enough to require corrections to the 
picture of isolated atoms, but not so much as to render the atomic descrip- 
tion completely irrelevant. This method, when not applied in oversimplified 
forms, describe electron propagation in any type of crystal (metal, semicon- 
ductors and insulators). Indeed, the tight binding model has a long history 
and has been applied in many ways and with many different purposes and 
different outcomes. Parts of the model can be filled in or extended by other 
kinds of calculations and models hke the nearly-free electron model. The 
model itself, or parts of it, can serve as the basis for other calculations. 
In the study of conductive polymers, organic semiconductors and molecular 
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electronics for example, tight binding like models are applied in which the 
role of the atoms in the original concept is replaced by the molecular orbitals 
of conjugated systems and where the inter atomic matrix elements are re- 
placed by inter or intra molecular hopping and tunneling parameters. 
For example, we imagine construction of a crystal from a hypothetical pe- 
riodic one- dimensional sequence of R equal atoms. In the case of negligible 
interaction among atoms, the same atomic orbitals centered in the different 
lattice sites would have the same energy; in the presence of interaction this 
R fold degeneracy is removed and evolves into an energy band. 
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Figure 3.1: Schematic representation of the crystal potential as a superposition 
of atomic-hke potentials, centered at the lattice sites tn- In the tight-binding 
approximation, the interaction between nearby atomic-hke orbitals of energy Ea 
leads to the formation of energy bands. Figure taken from [18J. 



3.2 The model 

In this section, we will introduce the tight-binding model that describes the 
system we have considered; we will introduce the general model ( the more 
general model that we will solve in this thesis) and the simplified model that 
will be solved in this chapter. 
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Figure 3.2: The star graph system (closed system). 



3.2.1 General model 



We consider h chains of sites of lengths Na, a = 1, . . . ,h.] the chains have a 
common vertex at the point O. In each chain, a particle can propagate with 
the hopping amplitude Va to closest neighbor^ The channels are coupled 
through the origin with hopping amplitudes taVa- We name this system 
Sh{Ni, . . . , Nh) system; indeed, this is a star graph. 
The total wave function of a particle in such a closed system is 



h Na 

J]J]C;^|a;n)+Co|0). 

a=l n=l 



(3.1) 



We assume that the numeration of the sites in each chain starts, = 1, at 
the site closest to the origin, while the site \a] No) will be later coupled to 
the external world. The level energy in each site equals 



Substituting (3.1) in the stationary Schrodinger equation with energy E, we 



^ For simphcity, fully justified by the localized nature of atomic orbitals, the hopping 
integrals involving second or further apart neighbours are assumed negligible. 
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obtain the following equations for the state (3.1): 



VaC2 + taVaCo = {E - (3.2) 

t;.(Q\i + CU) = {E- ea)C: n = 2, . . . , iV, (3.3) 

while the central amplitude Cq satisfies 

h 

Y.taVaC'^ = {E - eo)Co, (3.4) 

a=l 

where Eq is a local level at the center. The boundary conditions at the outer 
ends of a closed system are C'^^^^ = 0. The only coupling between the 



chains comes through eq. (3.4). Furthermore, when we will consider an open 
system, the last site of each chain will be coupled to the external world with 
coupling amplitudes . . . , y/jh- 



3.2.2 Simplified model 

In this chapter, we will consider the case of h chains all of identical length 
Na = N and identical hopping amplitudes Va = v among the h chains. We 
name this system Sh{N) system. In this case, we have hN + 1 sites, including 
the vertex. 

Taking into account all these considerations, the tight-binding Hamiltonian 
H with nearest neighbour interaction can be written as 

H = + He (3.5) 

where 



h N 

' (a; n\ 



a=l 71=1 

h N-l h N 

+ v\a] n) {a; n + 1| + w|a; n){a;n — 1\ 



(3.6) 



a=l 71=1 a=l n=2 



and 



H, = J2tav{\0){a;l\ + \a;l){0\) (3.7) 



a=l 
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It is clear from Eqs.(3.6) and (3.7) that H^c gives us the tight-binding model 



for h decoupled chains of sites (plus the origin) whereas He represents the 
coupling between the origin and the first site of each chain. 
Moreover, we set the level energy in each site (also the origin) equal to 
^0 = = 0. Thus, the Hamiltonian H becomes 



h N-1 



h N 



H = v\a]n){a\n + 1| + |a; n){a]n — 1\ 



h 



a=l n=2 



(3.8) 



+ J2tav{\0){a;l\ + |a;l)(0|) 



a=l 



Furthermore, starting from Eq.(3.1), we have that the total wave function of 



a particle in such a closed system is 

h N 



(3.9) 



a=l n=l 



Taking into account Eqs.(3.2), (3.3) and (3.4), the amplitudes of the station- 



ary state (3.9) with energy E satisfy the following equations, 
vC2 + tavCo = ECl 



while the central amplitude Cq satisfies 

h 

Y,tavCl = ECo, 

a=l 



n 



(3.10) 
(3.11) 

(3.12) 



3.3 Solution of the model for a closed system 



In this section, we will solve the simplified model presented in Paragraph (3.2.2 ) 
for a closed system. The equations that we have to solve are 

h 



J2tavC'^ = ECo 

a=l 

VC^ + taVCo = EC^ 

V + c:_,) = EC: 

vC'}f_i = EC^ 



a=l,. . . , h 
a=l,. . . , h] n=2, 
a=l,. . . , h 



(3.13) 



A^-l 
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One can easily notice that we have hN + 1 equations, as expected for hN + 1 
sites. 



The system (3.13) allows us to calculate numerically the eigenvalues of the 
system; it is enough, in fact, diagonalize the matrix that represents the sys- 
tem in the basis of the unknowns {Cq, C^}. 



A numerical solution of the system (3.13) gives us two eigenvalues outside of 
the normal Block hand, while the others hN — 1 are enclosed within the Bloch 
band. These two eigenvalues outside of the normal Bloch band are a distinc- 
tive feature of our star graph system; in the linear Id-chains considered so far 
in the literature, in fact, all the eigenvalues were confined within the Bloch 



band. In Fig.(3.3) we show the density of states (DOS) for the 5*4(100) system 
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t^-l 

o 

S-H 



Density of states for the S^(IOO) system 
v=l; tj=l, 1^=2, 1^-2, 1^=3 

yi \ \ \ ^ I \ \ \ I I \ \ \ I \ \ \ I I \ \ \ I ^ \ \ \ r 
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Figure 3.3: We plot the number of states as a function of the energy of the system 
for the 54(100) system with ti = 1,^2 = 2,t3 = 2,^4 = 3 and v = 1. We can see 
that the states are all within the band of Bloch, [-2, +2] in this case, except two, 
which are indicated with green circles. The quantities plotted are dimensionless. 
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4 chains each of 100 sites plus the origin - with ti = l,t2 = 2, ^3 = 2, 



and f = 1; it is clear from Fig.(3.3) that we have only two states outside 
of the normal Bloch band. This result is not due to the particular choice 
of parameters, but it is a general feature of the star graph system due to 
the coupling of the different chains with the origin, i.e. to the term He of 



Eq. (3.7). Indeed, if we set He = 0, all the states are confined within the 
Bloch band. 

A numerical study of the eigenvectors of the star graph system leads us to 
find that the two states corresponding to the eigenvalues outside the Bloch 
band are localized at the origin, where the h branches cross, while all the 
others are extended, Bloch-like states. 

In the following paragraphs, we will analytically calculate the eigenvalues of 



the states outside of the normal Bloch band (see Paragraph 3.3.1) and we 



will show that these states are localized at the origin (see Paragraph 3.3.2). 



3.3.1 Eigenvalues 



As we said before, we numerically find that the solution of the system (3.13) 
gives us two eigenvalues outside of the normal Bloch band, while the others 
hN — 1 are enclosed within the Bloch band. 

Eigenvalues outside of the normal Bloch band 

Now, we want to compute analytically the two eigenvalues outside of the 
normal Bloch band; we are particularly interested in those states because we 
will demonstrate that they are the only localized states of the system Sh{N). 
First of all, we want to have a linear system where the unknowns are a- 
independent, in order to have a (A^ + 1) x (A^ + 1) tridiagonal system which is 



easier to solve than the system (3.13). Starting from Eq.(3.13), we multiply 



and divide the left-hand side of the first equation by taV and we divide all 
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the other equations by taV; after that, we get 

h 



a=l 



taV 



taV 



Cn 



W+1 
taV 



+ v 



taV 



ECo 

E 
E 
E 



taV 
taV 



a=l,. . . ,h 

a=l,. . . ,h] n=2,. . . , 
a=l,. . . ,h 

(3.14) 
I9I c*^ 

It is easy to show that all ^ for i = 1, . . . ,N are a-independent if E and 
Co are different form zero. We are not interested in these cases because they 
do not lead us to states outside of the normal Bloch band. 
So, defining 



C,,= 



taV 



AT 



(3.16) 



and 



tot 



a=l 



(3.17) 



We first prove this statement for the particular case of = 3. In this case, we have 

h 



a=l 

vC^ + tavCa = EC1 
V + CI) = EC^ 
vC^ = ECl 



(3.15) 



0- 



Considering the fact that E, v and Co are a-independent, 



Replacing backward from the last equation, the second equation becomes 

9L ( ^ "^-E 

we find that also ^ is a-independent. In the general case, we find ^fj\i{E,v) — Cq 
where /at is a function dependent only on E and v for the case of N sites; again, the 
quantity ^ is a-independent. Please note that this proof is not valid if i? = or Co = 0. 



39 



3.3. SOLUTION OF THE MODEL FOR A CLOSED SYSTEM 



the system (3.14) becomes 



vC2 + Co - ECi 

vCn+l + vCn-1 - ECn 

fCjv_i — ECn 










n=2,. .., A^-l 



(3.18) 



This is the tridiagonal system (where the unknowns are a-independent) that 
we wanted to achieve and now we want to solve. This is an homogeneous 
system of + 1 hnear equations. 

Assigning a null value to each variable we get the null solution which is ob- 
viously not normalizable. 

Furthermore, we know from linear algebra that a square homogeneous lin- 
ear system has solutions (other than the trivial solution) if and only if the 
determinant of the matrix of the coefficients of the system is zero, namely. 



/ -E Ttotv^ 



det M = det 



1 





V 



-E 





V 

-E 










V 
V 












-E 

V 





-E 



\ 






V 



(3.19) 



This is a (A^ -|- 1) x (A^ -|- 1) determinant that describes the Sh{N) system 
and allows us to calculate the energy levels of the aforementioned system. 
The problem, now, is to find the determinant of the matrix M.] this matrix 
is tridiagonal, then we can calculate quite easily the determinant. 
Applying Eq.([R2|f1to Eq.(|3.19[), we get 



detM 



-E 



N-l 



-V 





-E 
1 




(3.20) 



11 



^See Appendix B. 
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Now, we have to calculate a power of a square matrix: 

A powerful method is to use the eigenvalue decomposition of the matrix B. 
It is straightforward to find that the eigenvectors of B are 

= ""I"' ^ = ""I"' j (3.22) 

where 

Now, let P be the matrix with these eigenvectors as its columns: 

P = ( ) (3.24) 

We can also compute the invers^P"^ of the matrix V 

Thus, the matrix B can be written as follows: 

B = VVV-^ (3.27) 
where "D is a diagonal matrix with the eigenvalues of B on the main diagonal: 

V=\ ^1 3.28 



4 



The inversion of 2 x 2 matrices can be done as follows: 

. ^ ( a hY^ I ( d -h\ 
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Computing the power of the matrix B: 

= (VV) [V-^V) V {P-^V) . . . {P-^V) VV-^ (3.29) 

Please note that the power of the matrix V is easy to calculate since it in- 
volves the powers of a diagonal matrix only. 



We are ready to replace these relations in Eq.(3.20 ); after some calculation 

we get 

det 7W = 

{E(-)f {E^ + EE^^) - T,W) - {E(+)f [E^ + EE^-) - T,^,v^) 



E(-) - E(+) ' ° 

(3.30) 

This is a polynomial equation that allows us, at least in principle, to find 
the eigenvalues outside of the normal Bloch band for any in fact, we just 
need to find the largest and the smallest root of this polynomial. 
From the computational point of view, however, is very difficult to provide 
an analytical solution with increasing degree of the polynomial to be solved; 
in fact, finding roots of certain functions, especially polynomials, frequently 
requires the use of standard numerical methods (such as Newton's method. 
Bisection method. Secant method, etc.). 

Here, we report the analytical solutions for some small N and for large N. 
Small N 

Where possible, we found analytically the roots of the polynomial of Equation 



(3.30); these are the results. 



For = 1: 



E = ±{^/Ttot)v (3.31) 



^We remind the reader that N is the number of sites for each chain, then the solution 
for N is the solution of the system that we named Sh{N) system. 
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For N = 2: 



E = ±{^/l + Ttot)v (3.32) 

For = 3: 



^2 + T,,, + V4T7i; 
E = ±\ ^ I V (3.33) 



For iV = 4: 



After these values of N, it is very difficult to find an analytical solution 
because of the complexity of the polynomial; however, is always possible to 
calculate a numerical solution with the above-mentioned numerical methods. 



Large N 



Here, we find the analytical solution for large N; this limit is very interesting 
because it appears very frequently in experimental situations. 



In this limit, from Eq. (3.30), after some calculations and remembering the 



definitions of E^ ^ and E^^\ we find that the equations to be solved are the 
following: 



E'^ + EVE^ - 4i;2 - 2Ttotv^ = 
E^ - EVE^ - 4i;2 - 2TtotV^ = 



(3.35) 
(3.36) 



or. 



E = ± 



tot 



VTt 



tot 



where Tm = ^ (^a) 



a=l 



(3.37) 



This relation gives the asymptotic eigenvalues for the system Sh as a function 
of Ttot and v; we remind the reader that v is the coupling between the sites 
of the same chain. In the limit of large Ttot, for the two states outside Bloch 
band we have E ^ ±\/T\^v, while all the others are always within the Bloch 
band [-2v,2v]. This is very important because one can control (by using Ttot 
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and thus also h) the energy distance between the states outside Bloch band 
and the others confined in this band. 



In Fig. (3.4) and Fig. (3.5) we show the results obtained with the 5*4(30) sys- 
tem; in these graphs, we compare numerical results (circles) with theoretical 



predictions (lines) given by the Equation (3.37). We plot only the positive 
eigenvalues outside of the normal Bloch band; the negative ones are symmet- 
ric with respect to the horizontal axis. 



In Fig. (3.4), we plot the eigenvalues as a function of v for different Ttof, we ob- 
tain three different straight lines with slope , according to Eq.(3.37). 
In Fig. (3.5), we plot the eigenvalues as a function of Ttot for two different 
values of f . In both figures, we can see the perfect agreement between the 
theoretical and numerical data; this is an evidence of the goodness of our 
analjd^ical model. Please also note that the theoretical data used for com- 
parison are those of the asymptotic regime; this means that for the 5*4(30) 
system we can consider the asymptotic regime already reached. 
We conclude stressing the fact that all the eigenvalues that we found here 
(for both small and large N) are a distinctive feature of our crossing config- 
uration system because they are outside of the normal Bloch band; in the 
linear Id-chains considered so far in the literature, in fact, all the eigenvalues 
were confined within the Bloch band. 



3.3.2 Eigenvectors 

A numerical study of the eigenvectors of the star graph system tells us that 
there are two types of eigenvectors: the two eigenvectors corresponding to 
the two eigenvalues outside of the normal Bloch band are localized at the 
origin, where the branches cross, while the other are extended, Bloch-like 
states. 

In this section, we will prove analitically that the two states outside of the 
normal Bloch band are localized and we will provide an analytical expression 
of their localization lengths in the asymptotic case. 



44 



3.3. SOLUTION OF THE MODEL FOR A CLOSED SYSTEM 



Eigenvalues as a function of the coupling between the sites of the same chain 

S_^{30) system 




V 



Figure 3.4: We plot one of the two eigenvalues outside of the normal Bloch band 
of the 5*4 (30) system as a function of v with different ta- In particular: 
Data set (a): ti = 15, t2 = 3, = 2, =1 =^ Ttot = 239 
Data set (b): ti = 7, = 5, ts = 4, t4 = 2 ^ Ttot = 93 
Data set (c): ti = 2, t2 = ^3 = 1) ^4 = 1 =^ Ttot = 7 

The circles refer to numerical data, while the lines represent theoretical data ob- 



tained by using Eq.(3.37) for different values of Ttot- The quantities plotted are 
dimensionless. 



Eigenvectors of the states outside of the normal Bloch band 

Now, we want to calculate the eigenvectors outside of the normal Bloch band. 



The linear system that we have to solve is Eq.(3.18) which is given here for 
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Eigenvalues as a function of the coupling with the origin 
84(30) system 




^2 10 20 30 ' 40 ' 50 ' 60 70 80 90 lOo' 1 lo 120 



T 



tot 



Figure 3.5: We plot one of the two eigenvalues outside of the normal Bloch band of 
the £'4(30) system as a function of Ttot for two systems with different v {v = 1 and 
V = 2). The circles refer to numerical data, while the lines represent theoretical 
data obtained by using Eq.(3.37) with v = 1 and v = 2. The quantities plotted 
are dimensionless. 



convenience: 



where 



vCr, 



TtotV^C^ 


= ECo 


VC2 + Co 


= ECi 


1 + vCn-l 


= ECn 


vCn-1 


= ECn 


CI 




taV 





n=2,. . . , N-1 



AT 



(3.38) 



(3.39) 
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anc 



E = ± 



tot 



(3.40) 



For convenience, which will be clear later, we make the following transfor- 
mation: 



a 



c 



N- 



with i = N 



(3.41) 



Thus, the system (3.38) becomes 



vCi = ECo 
vCo + vC2 = ECi 
vCi + vCs 



vCn-z + vCn-i 



EC2 

ECn^2 

ECn^i 



(3.42) 



Ttotv Cn-1 — ECn 
From the second row of the above equation, we obtain that 

C2 

and from the third row 



E ~ 

— Ci — Co 

V 



C, 



— C2 — C*i. 

V 



(3.43) 



(3.44) 



It is clear that the process can be iterated, and this leads us to recognize the 
following recurrence relatior^ 



E ~ 

Ck = — Ck-i — Ck-2 

V 



with k = 2 N -1 



(3.45) 



^Here, we are considering the asymptotic regime and then we use the asymptotic ex- 
pression for the eigenvalues; however, if one wants to consider a case with smaU N, one 
can get the result for that system simply using the appropriate eigenvalue calculated from 



Eq.(3.30l instead of Eq.(3.40|. 

For a brief overview of recurrence relations see Appendix A. In addition, please note 



that the recurrence relation (3.451 is valid until the third from last equation in (3.42 1 
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This is a linear homogeneous recurrence relation with constant coefficients 
of order 2. 

First of all, we determine the characteristic polynomial p{t) of this recurrence 
relation: 

' ' (3.46) 



p{t) = r 1 + 1 

V 

The roots of this polynomial are two, each of multiplicity 1: 

E ± - 4i;2 



£4 



2v 



(3.47) 



Therefore, we can write the general solution to our recurrence relation: 



with k = - 1 



(3.48) 



which gives us a closed-form expression for Ck- This is the most general 
solution; the two constants A and B can be chosen based on two given initial 
conditions Cq and Ci to produce a specific solution. 
We choose the following initial conditions for our recurrence relation: 

E 



Co 



(3.49) 



Obviously, these conditions do not return normalized eigenvectors, but we are 
interested in the ratios between the eigenvectors and those with conditions 



(3.49) are preserved. 

Thus, the system that we have to solve to find A and B is 

Co = A + B = l 



C, = Ai£_) + Bi£+) 



E 



(3.50) 



that gives, 
A 



E 

V 



£, 



B 



E 

V 



£- 



(3.51) 



Putting Equation (3.51) into (3.48) and after some algebraic manipulation, 

we get 



Ck oc 



VE^ - 4t;2 I 



{£- 



{£, 



k = A^- 1 



(3.52) 
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Recalling the transformation (3.41), we have 

^N-k OC — ===== 

VE^ - 4f 2 
or equivalently 



- {£+)"+' 



k = 0,...,N -1 (3.53) 



Ci OC 



vi-1) 



N-i 



(^4 



i = l,...,N (3.54) 



- 4t;2 

Now, we can compute the square of the absolute value of the ratio between 
the eigenvectors 



i = l,...,N -I (3.55) 



Up to now, we have not yet used the fact that we are considering h chains 
that are coupled by the following equation: 





2 




2 


{£^f-' - {£+f-' 













where 



E = ± 



tot 



tot 



(3.56) 



(3.57) 



One can check, by using Equation (3.55), taking into account (3.57) and 



the second row of (3.18), that the Equation (3.56) is automatically verified; 
therefore, we can say that the information on the coupling between the chains 
is contained only in the eigenvalues of the system. In addition, we note that 



if we substitute the Bloch eigenvalues into Eq.(3.55), we get the proper ratios 
between the Bloch eigenvectors; this is due to the fact that Eq.(3.55) is valid 
also for chains not coupled to the origin. 



Coming back to Equation (3.55), we want to compute explicitly £± for each 



eigenvalue; so, we must distinguish two cases, one for each eigenvalue. 
Considering E = ( = ) v and using the definitions of £±, Eq.(|3.47|) 



we get 



VTt 



tot 



£+ = V^tot — 1 £- 



VTt 



with Ttnt > 2 



(3.58) 



tot 
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On the other hand, with E = — 
1 



tot 



VTt 



tot 



V, we find 



tot 



1 



£- = VTtot 



with Ttnt > 2 



(3.59) 



Hence, substituting Eqs.(3.58) and (3.59) into Eq.(3.55), we obtain, regard- 
less of the considered eigenvalue, the following relation: 

2 / 1 \ 

B{Zot,N,t) 



tot 



iV- 1 



(3.60) 



where we have defined 



B{Ttot,N,i) 



1-(T, 



tot 



1 



,-{N-i) 



(3.61) 



1 - {T,^t - l)-(^+^-) 

This means that the ratio between two subsequent eigenvectors depends on 
z; however, this dependence is confined to the term B which is a boundary 
effect different from 1 only for i very close to A^. Neglecting the boundary 



effect^ namely setting B = 1, from Eq.(3.60) we get: 




1, 



.N 



(3.62) 



where i is the distance from the origin, where the h branches cross and the 
localization length ^ is equal to 



e = HTtot - 1) 



Ttot > 2 



(3.63) 



Therefore, the result indicates that the two states outside of the normal Bloch 



hand are localized with localization length ^, see Fig.(3.6) and Fig.(3.7). This 
result is a distinctive feature of our crossing configuration system. From 



Eq.(3.63), it is clear that the higher is Ttot, the higher is ^, which means that 
the two localized states are more localized. Please also note that this result 
is v-independent. 

It is also interesting compute the ratio between Co and C": 



^However, these effects will be crucial in determining the lifetime of localized states, 
see Chapter |4j Eq.( [4l8| ). 
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Localization length for systems with different different coupling with the origin 



T, = 33 

lot 




U 



T - 115 

tot 




Numerical data 
— Theoretical data 



Figure 3.6: In the simulation, we used a finite system consisting of 121 sites 



54(30) system. Theoretical data were calcolated using Equation (3.66). One can 



see the perfect correspondence between numerical results and theoretical predic- 
tions. We plotted only the first 20 value because the others were lower that the 
machine zero. The quantities plotted are dimensionless. 



ti 



Ttot > 2 



(3.64) 



tot 



This relation is, in general, different from one chain to another because it 
depends on the coupling ta with the origin of the chain considered; the higher 
is the coupling with the origin, the higher is this ratio, which means that the 
probability of being in that chain is higher. Thus, since the localization length 
is independent of the chain, the Cf give us h different sequences, where, 
however, the ratio of two successive terms is the same in all h sequences (i.e. 
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Localization length as a function of the coupling with the origin 



S^(30) system 




Figure 3.7: We show the locahzation length of the 54(30) system as a function of 
Ttot- The circles refer to numerical data, while the line represents theoretical data 



obtained by using Equation (3.66). The quantities plotted are dimensionless. 
chains). Namely, 





2 


Cf 


CI 


2 


C^ 


2 




Co 




CI 


C~o 




Co 







that gives us 



Co 



where 



and 



e = HTt,t - 1) 



tot) 



Ttot - 1 



(3.65) 



with i = l,...,N; Ttot > 2 (3.66) 



Ttot > 2 



Ttot > 2 



(3.67) 



(3.68) 
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3.4 An explanatory example 



In this Section, we present an example, in order to illustrate in a particular 
case what has been done in the previously. Let us consider h chains coupled 
with the origin with the same hopping amplitudes v. 
Namely, 



ta = 1 



with a = 1, . . . ,h 



(3.69) 



For example, a physical system that can be represented by this example is 



the crossing of h identical quantum wires. From Eq.(3.69), we have 

h 



tot 



EC. 

a=l 



h 



(3.70) 



In this paragraph, we will find the eigenvalues and the eigenvectors in the 
asymptotic regime. 



Eigenvalues - asymptotic regime 

A numerical solution of the system gives us two eigenvalues outside of the 
normal Bloch band, which is [— 2f , 2f], while the others are within the Bloch 



band; the density of states in a particular case is shown in Fig. (3. 8). Now 



we will find an analitical expression for these two eiegnvalues, starting from 
what has been done in the previous section. 



Substituting Eq.(3.70) into Eq.(3.37) we get 



E 



hv 



Vh 



with h > 2 



(3.71) 



The value h = 1 is meaningless; in fact, this would mean that there is a 
branch to the left of the origin, while on the right there is nothing. 
So, the first value that h can assume is h = 2; this means that we have two 
branches that cross in the origin, so we have a 1-d chain. For this particular 
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Density of states for the S jq(200) system 



100 



[/I 

a 
•t—i 



I 



10 



State outside 
Bloch band 




-1 
Energy 




Figure 3.8: We plot the number of states as a function of the energy of the system 
for the S'io(200) system with = 1 Va and v = 1. We can see that the states are 
all within the band of Bloch, [-2, +2] in this case, except two, which are indicated 
with green circles. The quantities plotted are dimensionless. 

case, we knowj^from the literature that the largest and the smallest eigenvalue 

^Consider a 1-d chain with nearest neighbour interactions composed by a finite number 
of sites from n — 1, . . . ,N with zero boundary conditions at the ends, n = and n = 1, 
where the chain terminates. The dispersion relation take the form 

nq 



E{q) 



Sq + 2v COS ^ ^ 



1 



1,2,. ..,N 



(3.72) 



In our treatment, we set — 0. Furthermore, we are interested in an infinite chain. It 
can be easily shown that for such a system the dispersion relation is 



E{fJ.) — 2v cos (ka) 



where 



TT , TT 

- <k< - 
a a 



(3.73) 



where a is the lattice spacing and fc is a continuous quantum number. From Eq.( 3.73"| , it 
is immediate to realize that, in the case of an infinite chain, the largest and smallest eigen- 
value are +2v and —2v, respectively. To ensure the convergence to the desired continuum 
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of a 1-d chain are respectively +2v and — 2f ; this is the same result that we 



have from Equation (3.71) setting h = 2. In the limit of large h, we have, for 
the two localized states, E ^ ±\/hv, while all the others are always confined 
in the Bloch band [-2v,2v]. Then, one can control, in a straightforward way, 
the energy distance between the two localized states and all the others by 
using the parameter h, i.e. the number of branches. 



Eigenvectors - asymptotic regime 

A numerical solution of the system leads us to find two classes of eigenvectors: 
the two eigenvectors corresponding to eigenvalues outside Bloch band are 
localized in the origin while all the others are extended, Bloch-like states. 
Taking into account what have been done in the previous section, we will 



compute the localization lentghs of the two localized states. From Eqs.(3.66) 
(3.67) and (3.68), considering ta = I Va, we get: 







2 


a 


2 




Co 




Co 





with i = 0, . . . , N] a = 1, 



h 



where 



e = ln(/i-l) 



h>2 



(3.74) 



(3.75) 



and where we have defined = Ci Va according to the symmetry of the 
problem. Therefore, we have h identical sequences, a-independent. 



In figure (3.9), we plot the localization length of 5/1(30) systems as a function 
of the number of branches; we set ta = 1 Va and v = 1. The line represents 



theoretical data obtained using Eq.(3.75) while the circles stand for numer- 



ical data; we can see the perfect agreement between the two series of data. 
Please note that theoretical data refer to the asymptotic regime; this implies 
that the system in question has already been achieved, with excellent ap- 
proximation, this regime. 

limit of a quantum wire, one has to let the lattice spacing a go to zero while increasing 
N such that their product (the length of the sample) remains constant. For more details, 
see Ref.[?]. 
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Localization length as a function of the number of branches 

S|^(30) system 




Figure 3.9: We plot the locahzation length of ^/^(SO) system as a function of h; 
we set = 1 Va and v = 1. The line represents theoretical predictions obtained 



using Eq.(3.75) while the circles represent numerical data; we can see the perfect 
agreement between numerical and theoretical data. The quantities plotted are 
dimensionless. 



3.5 An alternative derivation of the asymptotic 
regime 

In this section we will provide an alternative solution for the Sh system, which 
is valid only for the asymptotic regime and where we will neglect boundary 
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conditions. We start from Eq.(3.13) which is reported here for convenience: 

h 



Y,tavCl = ECo 



a=l 



vC^ + taVCo = ECl 



EC% 



a=l,. 
a=l 
a=l,. 



h 



h 

h; n=l, 



(3.76) 



Now, in order to have a system whose unknowns are a-independent, we define 

f 

and 

h 



(3.77) 



a=l 



tot 



(3.78) 



Thus, the system (3.76) becomes 






TtotV^C^ = 


ECo 




VC2 + vCq = 


ECi 


< 


vCn+1 + vCn-1 = 


ECn 




vCn-1 = 


ECn 



n=2. 



N-1 



(3.79) 



To solve this system we use the following "educated guess" , or an5at2p^ 



-kj 



with j = 0, 



N 



(3.80) 



fully justified by the fact that we are looking for states localized at the origin. 



where the branches cross. Substituting Eq.(3.80) into system (3.79), we get 

Zotve-' = E 

Ee-^ 



ve + v 



ve ^ + ve^ = E 



n=2,..., 



(3.81) 



ve 



E 



39 



^°Please note that the Ci are a-independent, see pag 

In physics and mathematics, an ansatz is an educated guess that is verified later 
by its results. An ansatz is the establishment of the starting equation(s) describing a 
mathematical or physical problem. 
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A solutionp^ of the system (3.81) leads us to: 

1 



and 



k 



E = ± 



HTt 



tot 



V. 



tot 



(3.82) 



(3.83) 



These are the same results for the asymptotic regime - and neglecting bound- 
ary conditions - calculated in the previous sections; see Eqs.(3.63)^^ and 



(3.37) respectively. 

Now a question arises: what's more in the approach presented in the previ- 
ous sections? First of all, the method described in previous sections allows 
us to calculate the energies of the localized states for any finite N by using 



Eq.(3.30) and not only for the asymptotic regime, as calculated in this sec- 
tion. Therefore, one can compute case by case the gap between the energy 
corresponding to the finite value of N and that given by the asymptotic limit, 
in order to understand when one can consider the asymptotic regime reached 
in relation to the accuracy that one wants to achieve. Secondly, it is also 
possible to calculate the localization length for any finite A^, by replacing the 



energy calculated with Eq.(3.30) in Eq.(3.55). Thirdly, the method described 



in previous sections takes into account boundary conditions, see Eqs.(3.60) 



and (3.61), while the one used in this section does not allows that; there- 
fore, the last method introduced considers only the case i3 = 1 in reference 



to Equation (3.60). We stress the fact that boundary conditions are crucial 
for the determination of lifetimes of localized states because the correction 
provided by B is very important at the last site that it is the site that must 
be considered for the lifetime, see Chap|4| Eq.(4.18). Finally, in the method 
described in previous sections we did not make any assumption on the so- 
lutions to be found, contrary to what we have done in the last section by 
choosing a suitable ansatz (not suitable if one takes into account boundary 
conditions). 

"'^^Solving the system we neglected boundary conditions. 
i^Provided to define £^ — 2k because of different definitions of k and 
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I'm personally convinced that computer science has a lot 
in common with physics. Both are about how the world 
works at a rather fundamental level. The difference, of 
course, is that while in physics you're supposed to figure 
out how the world is made up, in computer science you 
create the world. Within the confines of the computer, 
you're the creator. You get to ultimately control every- 
thing that happens. If you're good enough, you can be 
God. On a small scale. 



Linus Torvalds - creator of Linux 




Open tight-binding model: 
transition to superradiance 



4.1 Effective Hamiltonian for the star graph sys- 
tem 

Let us consider a system composed of h chains of sites of lengths Na, a = 
1, . . . ,h] the chains have a common vertex at the point O, the origin. We 
named this system Sh{Ni, . . . , Na) systen^ see Chapjsj Sec. (3. 2). The total 
wave function of a particle in such a closed system is 

h Na 

i^) = EE^«i«'^)+^o|o)- (4.1) 

a=l n=l 



^Furthermore, if we set Na ^ N we obtain the Sh{N) system. 
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V t^v vjy 




Figure 4.1: The star graph system. 



We assume that the numeration of the sites in each chain starts, n = 1, at 
the site closest to the origin, while the site \a;Na) will be coupled to the 
external world. The level energy in each site (also the origin) equals zero, 
for simplicity. Thus, the tight-binding Hamiltonian with nearest neighbour 
interaction of the Sh{Ni, . . . , A^^) system can be written as 



h Na-l h Na 

H = v\a; n) {a; n + 1| + f |a; n){a;n— 1| 

a=l n=l a=l n=2 

h 

+ I]M(|0)(a;l| + |a; 1)(0|) 

a=l 



(4.2) 



This is the Hamiltonian for the closed system; now, we want to introduce 
the coupling with the environment. The outside world is characterized by a 
continuum of states |c, E) where c = 1, . . . , M is a discrete number labeling 
M channels and E is a continuum quantum number representing the energy. 
In this thesis, we will couple the last site of each chain to a single external 
channel; this is a choice because it is possible to build more complicated 
couplings with the outside world, that leads to hierarchical models, see for 
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example [19]. Under these considerations, the channel states can be written 



as \a, E) where a = l,...,/i, according to Eq.(4.2). Taking into account 



Eq.(2.42), we have 



n = H-lw W^^ = J2KK* (4.3) 



2 

a=l 



where H is the same defined in Eq.(4.2) and A'^ is the transition amplitude 



between the intrinsic state \a;n) and the continuum \a,E). For system in- 
variant under time reversal, both H and W are real symmetric matrices, and 
the coupling amplitudes A'^ can be taken as real. Furthermore, we set A'^ all 
equal to ^/7; so, 7 is a parameter that controls the coupling strength with 
the external world. Due to the coupling to the channel states, the states 
\a,Na), a = 1, . . . acquire a finite width 7. Now, we can write the full 
effective Hamiltonian for this system: 

h 

^ = H-^7 5Z|a,iVa)(a,iVa| (4.4) 

a=l 



The diagonalization of the non-Hermitian effective Hamiltonian (4.4) gives 
us the complex eigenvalue^ 

£i = E,- (4.5) 

corresponding to resonances centered at Ei with widths Fj that determine 
the lifetime of a state, ~ h/Ti. 

Imaginary part of the eigenvalues outside Bloch band 

Now, we want to provide an estimation of the imaginary part Tioc of the 
eigenvalues Sioc outside Bloch band. 

First of all, we define the wave function for the localized states as 

h N 

l^ioc) = E E ^tioch n) + Co,ioc\0) (4.6) 

a=l 1=1 



These complex eigenvalues of H coincide with the poles of the S-matrix, see Chap§ 
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and defining Tfot = Sa^a' have for tlie asymptotic regim^ 

1 



tot 



B{Ttot,N,t) 



i = l,...,N -1 



witli 



B{Ttot,N,t) = 



1 - (Ttot - 1 



,-{N-t) 



(4.7) 



(4.8) 



1 - mot - 1)-(^+^-^) 

Tlie term i3 is a boundary effect that it is very important (namely, very 
different from 1) for i very close to which is precisely the case that must 
be considered in this paragraph. Let us now consider the ratio 

2 

(4.9) 



'^N,loc 


2 




2 


'^2,loc 


2 




2 


^N,loc 










^l,loc 




^2,loc 




^N-l,loc 



that, taking into account Eqs.(4.7), (4.8) and the fact thal|^ 



r^a 2 

Co,loc 

lead us to the following relation: 

t2 



2 

'^N,loc 



O,loc 



Ttot - 1 



N-l 



(4.10) 



tot 



i=i 



(4.11) 



which, using the definition of ^, ^ = ln{Ttot — 1), can be written as 



r^a 2 
'^NAoc 



Co, 



loc 



N-l 



(4.12) 



i=l 



Furthermore, we know from the literatur^ that, for small 7, the width of a 
state can be written as 



rioc = X^|(a;Ar,|V'Zoc)r7 



(4.13) 



a=l 



Eqs. 



3.601 and (3.611. 



^See Chap 

''See Chap^ Eq.(|^| 
^See for example Refe![20i and |2|]. 
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and it agrees with the intuitive argument that the decay executes the decom- 
position of the state, isolating the components matched to the specific open 
channels. For semplicity, we consider h chains of identical length; hence, 



Eq.(4.13) becomes 



Tioc = Y.\{a-N\i;ioc)\- 



'7 



(4.14) 



a=l 



Taking into account Eqs.(4.6) and (4.12), we obtain 



N-l h 
a=l 



1=1 



that gives us 



Af-l 



or, using the definition of ^, 



(4.15) 



(4.16) 



i=l 



N-1 



loc 



tot 



^) 1=1 



(4.17) 



Moreover, we know that the lifetime of a state is Tj ~ h/Tf, so, setting h = 1, 
we get for the localized states: 




(4.18) 



which is valid only for small 7. We can see from equation (4.18) that the 



larger (the number of sites for each branch), the greater the lifetime of 
the localized states; hence, for very large A^, we have tioc — for small 



values of the external coupling 7. Moreover, we note that in Eq.(4.18) the 
factor |Co,iocP appears and it can be calculated in a straightforward manner 
from the knowledge of the relationship between two successive eigenvectors. 



Eq.(4.7), and by imposing the normalization of the wave function. 



In Fig. (4.2), we show the imaginary part Tioc of one of the two eigenvalues 



corresponding to the localized states as a function of the external coupling 7 



63 



4.1. EFFECTIVE HAMILTONIAN FOR THE STAR GRAPH 

SYSTEM 



(the other is the same). We considered the 5*4(10) system setting = 1 Va 
and V = 1. The sohd black hne refers to theoretical data obtained by using 



Eq.(4.17) while circles refer to numerical data; for small 7, there is perfect 



correspondence between theoretical and numerical data. 



Imaginary part of one of the two eigenvalues outside Bloch band 
(localized states) 



° Numerical data 
— Theoretical data 




Figure 4.2: We plot the imaginary part Tioc of one of the two eigenvalues outside 
Bloch band corresponding to the localized states as a function of the external coupling 
7. We considered the 54(10) system with = 1 Va and v = 1. The solid line represents 



theoretical predictions obtained by using Eq.(4.17l while circles refer to numerical data; for 
small 7, we can see the agreement between numerical and theoretical data. The quantities 
plotted are dimcnsionless. 
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4.2 Superradiance transition 



4.2.1 Evolution of the complex eigenvalues 



Evolution of the eigenvalues of the system (5) 




y 

Figure 4.3: The evolution of the eigenvalues £i = Ei — z/2r,j of the effective Hamiltonian 
is shown as the external coupling 7 is varied. The 6*4(5) system (2f sites) with ta — l^a is 
considered. The quantities plotted are dimensionless. (a) The real part of the eigenvalues 
£i as a function of the coupling parameter 7. (b) The imaginary part of the eigenvalues £i 
as a function of the coupling parameter 7. The emergence of four superradiance states is 
clearly visible: one state with multiplicity 3 - red circles - and one state with multiplicity 
1 - blue circles; see text for further explanation. 

Now, we want to show how and when the superradiance transition occurs. 
First of all, we know that the eigenvalues of the effective Hamiltonian l-L can 
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be written as Si = Ei — i/2Ti. For 7^1, the real part of the complex 
eigenvalues, E^, are close to the eigenvalues of the internal Hamiltonian H 
and the imaginary part, Fj, gives the width of the isolated resonances. In the 



upper panel of Fig. (4.3), we show Ei as a function of the coupling constant 
7 for the 54(5), which stands for 4 chains of 5 sites each, with ta = 1 Va; 
similarly, in the lower panel, we plot F, as a function of the coupling 7. We 
see that, as 7 increases, h = 4 states - blue and red circles in the lower 



panel of Fig. (4. 3) - acquire the most of the width and become short-lived or 



superradiant. Still looking at the lower panel of Fig.(4.3), we see that there 
are two states that have the same value of F up to a critical value 7 ~ 3 
where a sort of bifurcation occurs and one state becomes superradiant at 
the cost of the other. At exactly the same critical value of 7, the real part 
of these two states both go to zero, as can be seen in the upper panel of 
the figure. The superradiant state will appear only at the critical coupling, 
7^3. The other three (degenerate) superradiant states (red circles in 
the lower panel) are corresponding to Ei = 0. Summarizing, we can state 
the following: as the coupling increases and reaches a critical value, the 
resonances overlap, and a sharp restructuring of the system occurs. Beyond 
this critical value, a few resonances become short-lived states, leaving all 
other (long-lived) states effectively decoupled from the environment. This 
general phenomenon is referred to as the superradiance transition, due to 
its analogy with Dicke superradiance in quantum optics, see Ref.[T5]. At 
the critical point a transition takes place which is caused by the feedback 
between environment and system. 



In Fig. (4.4), the evolution in the complex plane of the eigenvalues of the 
effective Hamiltonian is shown as the external coupling 7 is varied. Arrows 
indicate the motion of the poles in the complex plane as 7 is changed. The 
5*50(5) system (251 sites) with = 1 Va and t> = 1 is considered. According 
to Chapj3]-Sec.(3.3), for the closed system (7 = 0) we have two eigenvalues 
outside of a normal Bloch band while the others 249 are enclosed within 



the Bloch band. In the panel (a) of Fig. (4.4), we plot the evolution of the 
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Evolution of the eigenvalues of the system (5) in the complex plane 
as the external coupling y is varied 
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Figure 4.4: The evolution in the complex plane of the eigenvalues of the effective Hamil- 
tonian is shown as the external coupling 7 is varied. The >S'5o(5) system (251 sites) with 
ta — lya and u = 1 is considered. Arrows indicate the motion of the poles in the complex 
plane as the external coupling 7 is increased (0 < 7 < 10 in the simulation). The quantities 
plotted are dimensionless. (a) The evolution in the complex plane of the 249 eigenvalues 
within Bloch band. We reported on the symbol the multiplicity of the eigenvalue (blue). 
After a critical value of 7, the superradiance transition occurs. After that, we have 51 
states with Ei ~ 0; 50 states are superradiants and one state is not superradiant. (mul- 
tiplicity in red), (b) The evolution in the complex plane of one of the two eigenvalues 
outside Bloch band. 



eigenvalues inside Bloch band; we reported on the symbol the multiplicity of 
the eigenvalue (blue). Generally speaking, we can see that, as 7 increases, 
the poles tend to go toward the center of the band; moreover, it is interesting 
to look at what happens at the two states closest to E' = in more detail. 
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As 7 increases, these two states approach E = 0, until a critical value of 7, 
when the superradiance transition occurs. After that, we have 51 states with 
Ei ~ 0; 50 states are superradiants and one state is not superradiant. In 



the panel (a) of Fig. (4.4), we reported in red the multiplicity of the states 
with ~ after the superradiance transition. In the panel (b), we plot the 
eigenvalue trajectory of one of the two localized states outside Bloch band 
as a function of the overall coupling strength 7. We see that the eigenvalues 
outside the Bloch band in the complex plane are weakly dependent on the 
external coupling. 

4.2.2 Evaluation of critical 7 



As we explained in Chap|2| Sec. (2. 4), we can distinguish two different regime 



for increasing 7. At weak coupling, all internal states are similarly affected 
by the opening and acquire widths proportionally to 7. In the opposite limit 
of large 7, only h states (where h is the number of channels) will have a 
width proportional to 7, while the widths of the remaining states fall of 
as 1/7. In order to find the critical value of the parameter 7 at which the 
superradiance transition occurs, we may analyze the average width (F) of the 



N — h narrowest widths as a function of the coupling 7, see Fig.(4.5). At the 
critical value of 7, the average width (F) peaks and begins to decrease. We 
can evaluate this critical value of 7 using the criticality criterion discussed 



earlier in Chap|2| Eq.(2.57). Roughly, the transition occurs when (T) / D ^ 1 



where D is the mean level spacing of the Hamiltonian for the closed system. 
In our system, there is lot of degeneracy; for the Sh{N) system (i.e. hN + 1 
sites), we numerically find that there are M different energy levels, where 
M = 2N + 1; so, we define an effective mean level spacing Dejf'- 

^ 4v 4v 2v , 
Moreover, the average width is 
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Evaluation of critical y for different systems 

(a) (b) 

0.15R 




^0 1 2 3 4 5 6 

y/v 



1 2 3 4 5 6 

y/v 



Figure 4.5: The average width (F) is shown as a function of 7. When the number of 
sites is large, the transition to superradiance is shown to occur at ^cr/v — 2, regardless of 
the value of ft., see text. The solid line corresponds to an average over all K widths, while 
the symbols are obtained by averaging over the K — h smallest widths. On each panel, 
we indicated in the box the system under consideration; for all systems we set ta = 1 Va. 
The quantities plotted are dimensionless. 



Therefore, the criticality criterion becomes 



(r) 



7c 



(4.21) 



Deff 2v 

This imphes jcr = 2u; thus, the critical value of 7 is independent of h and 



ta Va. In Fig.(4.5), we show the average width (F) as a function of 7. The 
solid line corresponds to an average over all K widths, while the symbols 
are obtained by averaging over the K — h smallest widths. The numerical 
computation of the critical 7 indicates that the above theoretical estimation 
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is in perfect agreement with numerical data; thus, for the systems 5*4(40) 
and 5*10(40), we can consider the asymptotic regime aheady achieved, see 



panel (a) and (b) of Fig.(4.5). Furthermore, if we look at panel (c) and 
(d), we notice that the critical 7 is different from that of the asymptotic 
regime and also depends on the h value; these effects, however, are due 
to the finite range considered and disappear in the asymptotic regime. In 



Evaluation of critical y for the S^CTO) system for different value of t 




Figure 4.6: The average width (F) is shown as a function of 7. When the number of sites 
is large, the transition to superradiance is shown to occur at ^cr/v ~ 2, regardless of the 
value of ta- In the simulation we considered the 5*4(70) system with different values of tj,: 
ta = \ Va(green line) and ti = 200, t2 = 100, = 20, ^4 = 50 (red circles). We can see the 
perfect correspondence between these two sets of data. The solid black line corresponds 
to an average over all K widths, while the green line and the symbols are obtained by 
averaging over the K — A smallest widths. The quantities plotted are dimensionless. 

order to show that the critical 7, in the asymptotic regime, is independent 
of the choice of ta, we consider the 5*4(70) system with two different sets 
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of t„: in the first case, we choose t„ = 1 Va, while in the second we set 



ti = 200; t2 = 100; tg = 20; ^4 = 50. The result is shown in Fig.(|4^. The 
solid black line corresponds to an average over all K widths, while the green 
line and the symbols are obtained by averaging over the K — 4 smallest 
widths; in particular, the green line refers to the first case, ta = 1 Va, while 
the symbols refer to the second case. The perfect agreement between these 
two sets of data is evident. 

Chains with different lengths 

Now, we want to investigate the behavior of the critical 7 for a system whose 
chains have a different number of sites, in general A^^ 7^ Na' ii a ^ a' . This 
question is relevant for experiments where the exact number of sites can be 
hardly controlled; in fact, the sites of our model can represent, for example, 
biological molecules or ions and it is therefore very common have to deal with 
chains of different lengths. First of all, we study the asymptotic regime, i.e. 
when the total number of sites X]a^« large; after that, we will take a 

look at systems with a small number of sites. 

In order to shed light on asymptotic situation, we consider two different sys- 
tems: the 5*4(60) system and the 5*4(35, 70, 53, 82) system, both with = 1 
Va. These two systems have the same number of sites (241). In the panel 



(a) of Fig. (4.7), we show the average width (F) as a function of 7. The solid 
black line corresponds to an average over all K widths, the solid magenta line 
is obtained by averaging over the K — A smallest widths for the 5*4(60) system 
and the blue circles are found by averaging over the K — A smallest widths 
for the 54(35,70,53,82) system. We can see the perfect correspondence be- 
tween these two series of data: this is not due to the particular choice of the 
parameters Na, but is a general behavior of our system in the asymptotic 
regime. Thus, we can state that the critical 7, in the asymptotic regime, is 
independent of the values of Na- 

Now, let us have a look at what happens when the system is composed of 
a small number of sites. We consider three different systems each consisting 
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Evaluation of critical y for the S system for different value of N 
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Figure 4.7: The average width (F) is shown as a function of 7. When the total number 
of sites (X^a '^^) large, the superradiance transition is shown to occur at ^cr/v — 2, 
according to text, regardless of the different values of Na- In the simulation we considered 
the S4{Ni, N2, N3, N4) system with different values of Na and ta = I Va. The quantities 
plotted are dimensionless. (a) Na — 60 Va(solid magenta line) and A^i — 35, N2 — 70, 
N3 — 53, N4 = 82 (blue circles). We can see the perfect agreement between these two 
series of data. The solid black line corresponds to an average over all K widths, while the 
green line and the symbols are obtained by averaging over the K — A smallest widths, (b) 
iVa = 5 Va (red circles), iVi = 2, N2 = 6, N3 = 4, N4 ^ 82 (green circles) and A^i = 2, 
N2 — 2, N3 = 2, N4 = 14 (blue circles). The solid black line corresponds to an average 
over all K widths, while the symbols are obtained by averaging over the K — A smallest 
widths. We can see that, if the total number of sites is small, systems with different Na 
exhibit very different behaviors. 
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of 21 sites: ^4(5), ^4(2, 6, 4, 8) and ^4(2,2,2,14) all with t„ = 1 Va. In the 
panel (b) of Fig.( |4.7[ ), we plot the average width (F) as a function of 7. The 
solid black line corresponds to an average over all K widths, while symbols 
is obtained by averaging over the K — 4 smallest widths. It is evident that 
the average over the K — 4 smallest widths for the 5*4(5) system (red circles) 
shows a single peak: this is the unique critical value of 7. On the other hand, 
if we look at the two systems with different A^^^ (blue and green circles), we 
can easily see that the average over the K — A smallest widths exhibits more 
than one peak, so we have more than one critical 7; this feature is very evi- 
dent for the data series related to the 5*4(2,6,4,8) system (green circles). 
However, these effects are due to the small number of sites of the systems 
considered and disappear in the asymptotic hmit. 



4.2.3 Evolution of the localization length 



In Chapter Sec. (3.3), we showed that for the Sh{N) system we have two 



localized states corresponding to the two eigenvalues outside of a normal 
Bloch band. Writing the wave function for the localized states as 

h N 

li^ioc) = E E ^hoch n) + Co,ioc\0) (4.22) 

a=l n=l 

neglecting boundary effects (namely, setting B = 1) and defining T^t = 
Sa^a' have for the asymptotic regime: 



C. 



n.loc 



Co,loc 

where 



^(t„,Ttot)e-«("-^) withn = l,...,iV; T^^ > 2 (4.23) 



e = HTtot - 1) 



Ttnf. > 2 



(4.24) 



and 



Ttnt > 2 



(4.25) 



tot 
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Localization length as a function of y for the S (100) system 




10 20 30 40 50 
i 



^ (c) 

(N 0^-^^^ I 

U°-20 
B-40 
J -60, 



7=2 



10 20 30 40 50 

i r 




10 20 30 40 50 



(d) 



U°-20 
^^40 



-60 



Y=5 



• Numerical data (open system) 
— Theoretical data (closed system) 
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Figure 4.8: In the simulation, the 54(100) system with ia = 1 Va and i; = 1. Theoretical 



data (for a closed system) were calculated using equations (4.231, (4.24) and (4.24). One 
can see that the localization length is j -independent. 



Now, we want to study the behavior of the localization length as a function 



of the openness of the system. In Fig(4.8), we plot the evolution of the 



localization length of one of the two localized states as the external coupling 
7 is varied. The circles refer to numerical data, while the straight line refers 
to theoretical data (closed system); theoretical data were calculated using 
equations (4.23), (4.24) and ( 4.24[ ). We notice that the localization length 
is -independent. In the simulation we considered the 54(100) system with 
= 1 Va and v = 1. 

However, this result is general; it is also valid for different values of h, ta 



and V. Actually, from Fig. (4.4) on pag,67, we can see that the imaginary 



part of the eigenvalues outside Bloch band are very small, which means that 
localized states are very little affected by the opening of the system; in this 
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sense, the results shown in Fig. (4. 8) should not be considered as a surprise. 



4.3 Resonance structure 



In order to show the consequences of the superradiance transition on the 
transport properties, here we analyze the resonance structure, by considering 
the transmission T (E). In Chapj2| Eq.(2.49) and Eq.(2.53), we found that 
the transmission for the process from channel a to 6 is 



T'^^E) = IT"" (E) 



(4.26) 



where is the number of the internal states and Si = Ei — i/2Ti are the 
eigenvalues of the effective Hamiltonian which coincide with the poles of 
T'^'^^E). Now, we want to compute the trasmission from channel a to 6 for 
the Sh{N) system at the energies corresponding to localized states, Sioc = 
Eioc — if^Tioc] the last site of each chain is coupled to the external world with 
coupling amplitudes = y/^, Va = 1, . . . , /i. 
We define the wave functioij^ for the locahzed states as follows: 



N 



IKc) = EQo>5^)+^0,Zoc|0) 



(4.28) 



1=1 



and we compute the transmission at the energies corresponding to localized 
states: 



T"*(E = E, 



loc) 



{a;N\ij^J^{^t\b;N)^ 



Furthermore, from Eq.(4.12) we know that 



5-F 

)-■- loc 



N-1 



(4.29) 



(4.30) 



1=1 



^Please note that the wave function defined in (4.28 1 is not normalized; the normalized 
wave function is the following: 



h N 

\i^ioc)^J2J2^hoc\a;n)+Co,ioc\0) 

a=l i=l 

but here we consider IV'foc) convenience. 



(4.27) 
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where 



tot 



a=l 



e = ln(Ti„,-l) 



(4.31) 



and from this chapter, see Sec. (4.1) - Eq.(4.14), that 



(4.32) 



a=l 



Taking into account all these considerations and after some calculation, we 
get 

h 



T"''(E = E, 



loc ) 



IT |2 
I tot I 



tot 



a=l 



(4.33) 



This formula gives us the transmission from channel a to 6 at the energies 
corresponding to localized states. If we set, for example, ta = 1 Va, we have 
■jab _ /^^f^2 \Jq^^]j-^ ^]2g transmission is the same in all channels, as one 
might expect given the symmetry of the system]^ 

Numerical simulations show that the transmission calculated at the energies 
corresponding to localized states is the same for all the other states; see for 



example Fig. (4.9). 



Example 1: identical coupling between the chains and the origin 

In this paragraph, we will consider the particular case of ta = 1 Va; this 



means that all the chains are equally coupled with the origin. In Fig. (4. 9), 
the transmission T^^ is shown as a function of the energy E for several values 
of the external coupling 7 for the 6*4(5) system with = 1 Va and f = 1; 



from Eq.(4.33) we obtain that the transmission peaks reach the value of 1/4, 



see Fig.(4.9). Moreover, it is clear from Fig.(4.9) that the superradiance 



transition has a clear signature in the resonance structure. 

For the Sh{N) system (i.e. hN + 1 sites), we numerically find that there 

are 2N + 1 different energy levels because of degeneration; so, for the 5*4(5) 



^See ReLM- 
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Figure 4.9: The transmission T is shown as a function of energy E for the 5*4(5) system 
(21 sites) for different values of the external coupling 7. T'^'' is the transmission from 
channel a to b. In this figure, we plotted only the results for a — 1 and 6 = 2. The results 
for a — 2,3,4 and b = 1, 2, 3, 4 are the same. We considered v = 1. The quantities plotted 
are dimensionless. 

system, we have 11 different energy levels for the closed system. Hence, the 



resonances in Fig. (4.9 ) are 11: the two external narrow resonances correspond 
to the localized states, outside the normal Bloch band, while the others nine 
correspond to extended, Bloch-like states. 



For weak coupling, 7 = 0.3, see panel (a) of Fig.(4.9), there are 11 isolated 
resonances which corrispond to the 11 classes of different eigenvalues that ac- 
quire a small imaginary part. Since 7 is small, the imaginary part gives the 
width of the isolated resonances. The width decreases away from the band 
center {E = 0), reaching a minimum with the last two resonance peaks. 
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corresponding to the localized states; indeed, these resonance peaks are the 
narrowest. As we increase the external couphng, the transmission increases 
and near 7 = 1.5 all resonances (except the two outer corresponding to lo- 
calized states) start to overlap, see panel (b). 

At 7 = 2, see panel (c), the five central resonances are almost completely 
superimposed. For 7 = 5 (strong coupling), we notice, see panel (d), that 
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Figure 4.10: The evolution in the complex plane of the eigenvalues of the effective 
Hamiltonian is shown as the external coupling 7 is varied. The >5'4(5) system is considered. 
We reported on the symbol the multiplicity of the eigenvalue. The emergence of the 
four superradiance states is clearly visible in the panel (d). The quantities plotted are 
dimensionless. 

we return to the case of isolated resonances but this time, their number is 
reduced to 11 — 2 = 9 resonances; therefore we lose only two resonance peaks 
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4.3. RESONANCE STRUCTURE 



because of the degeneration of the system. Indeed, the peak at ~ is 
associated with three states and each of the two states closest to ~ has 
multiphcity 1; the superradiant states are four, two at ^ and the other 
two closer to -E ~ 0, as expected for four channels. Thus, the peak in the 
middle of the band remains because it had (before the superradiance transi- 



tion) multiplicity 3, while the other two neighbors disappear. In Fig. (4. 10), 
we show the evolution in the complex plane of the eigenvalues of this system 
as the external coupling 7 is varied. 



Furthermore, we can notice from panel (d) of Fig. (4.9) that the width of the 



resonances is smaller for energies closer to the band edge, just as it was for 



weak coupling. Finally, we note that the resonance peaks in the figure (4.9) 
reach the same value of 1/4 in each branch, as expected since the ta are all 
the same. 

Example 2: different coupling between the chains and the origin 

Let us consider an example where the branches have different couplings with 
the origin; here, we consider the ^'4(4) system setting ti = 3; ^2 = 2; = 
1; ^4 = 1 and v = 1. We plot the transmission from several input channels to 



different output channels for the states within Bloch band in Fig.(4.11) and 



for one of the two localized states in Fig. (4. 12); one can easily check that the 



resonance peaks reach the values found in Eq.(4.33) for all states. Looking at 



the figure (4.11 ), we see that after the superradiance transition two resonance 
peaks disappear; actually, the superradiant states are four because there are 
four open channels. Only two peaks are missing because of degeneration 



of our system; in other words, the situation is very similar to Fig.(4.10) 



related to the previous example. Finally, in Fig. (4. 12), we plot one of the 
two eigenvalues outside Bloch band; we notice that the resonances related to 
the localized states are very narrow compared with those corresponding to 
extended states and we can also appreciate the energy shift due to the fact 
that we considered three different external couplings, 7 = 0.3 (blue line), 
7 = 2 (green line), 7 = 20 (red line). 
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4.4. INTEGRATED TRANSMISSION 



Transmission for the S^(4) system for two differents y 

t, =3;!^ =2;^ = 1; 1^=1 
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Figure 4.11: The transmission T is shown as a function of energy E for the £'4(4) system 
(17 sites) for two different values of the external coupling 7, 7 = 0.3 (solid magenta line), 
7 = 5 (dashed blue line) before and after the superradiance transition, respectively. We 
set: ti = 3;t2 = 2; is = l;t4 = 1 and v = 1. Here, we plot only the states within Bloch 



band. T"'' is the transmission from channel a to b. It is easy to check that Eq.(4.33) is 
valid also for the states within Bloch band, [-2, +2] in this case. The quantities plotted are 
dimensionless. 



4.4 Integrated transmission 



It is also interesting compute the normahzed integrated transmission: 



AE 



T (E) dE 



AE E^nax EfYiin 



(4.34) 
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4.4. INTEGRATED TRANSMISSION 



Transmission for the S^(4) system - one of the two localized states 
t,=3;t^ =2;t3 = l;t^=l 




4.0085 4.009 
E 



4.0095 



Figure 4.12: The transmission T is shown as a function of energy E for the 5*4(4) 
system (f 7 sites) for different values of the external coupling 7 only for one of the two 
states outside Bloch band. T"^^ is the transmission from channel a to b. We set: ti — 
3; t2 = 2; = 1; <4 = 1 and w = 1. In this figure, we plotted only the results for a — 1 and 
b — 1. The quantities plotted are dimensionless. 



where the interval AE includes all the system eigenvalues. 



In Fig. (4. 13), the integrated transmission / is shown as a function of external 
coupling 7 for two different systems: 5*4(10) and 5*4(50) systems. We plot 
also the normalized integrated transmission for a one dimensional chain for 
comparison: 

X = / T{E)dE = , (4.35) 

"^vj ^ ^ 4 + (7/t;)' ^ ^ 

divided by 4. Indeed, since the system under consideration consists of four 

chains coupled in the origin, it is reasonable to compare the integrated trans- 



^See Ref.[T] for more details. 
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4.4. INTEGRATED TRANSMISSION 



Integrated transmission for different systems 




123456789 10 

y/v 



Figure 4.13: The integrated transmission / is shown as a function of external couphng 
7 for two different systems: 54(10) and 54(50) systems. 5"'' is the transmission from 
channel a to b. In this figure, we plotted only the results for a = 1 and b — 2. The 
results for a — 2, 3, A and b — 1, 2, 3, 4 are the same. We considered v = 1. The quantities 
plotted are dimensionless. The integrated transmission I^^ has a maximum at jtr/v ~ 2: 
this is the critical value of 7 for the transmission. This value is almost independent from 
the number of sites considered as we can see from the figure above. We plot also the 
theoretical integrated transmission for a chain of sites, divided by 4. 

mission with the results already known for the system consisting of one chain, 
divided by four. 

The integrated transmission J^^ has a maximum at '~ftr/v — 2: this is the crit- 
ical value of 7 for the transmission. This value is almost independent from 



the number of sites considered as we can see from figure (4.13). Moreover 



for large N, the transmission maximum is reached precisely at the superra- 



diance transition: — 'Jtr = 2f, compare Fig. (4.5) and Fig. (4. 13); so, the 
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4.4. INTEGRATED TRANSMISSION 



maximum of the integrated transmission is strongly related to the superra- 



diance transition. Moreover, we can see from Fig. (4. 13) that the normalized 
integrated transmission for a one dimensional chain of sites is very differ- 
ent from the normalized integrated transmission obtained with our crossing 
configuration system. 
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Fortunately science, like that nature to 
which it belongs, is neither limited by time 
nor by space. It belongs to the world, and 
is of no country and of no age. The more 
we know, the more we feel our ignorance; 
the more we feel how much remains un- 
known; and in philosophy, the sentiment of 
the Macedonian hero can never apply, there 
are always new worlds to conquer. 

Sir Humphry Davy - British chemist 




Perspective 



5.1 Disorder 



We know from previous chapters, see Eqs.(3.5), Eq.(3.6) and Eq.(3.7), that 



the tight-binding Hamiltonian H with nearest neighbour interaction for the 
star graph system Sh{N) can be written as 



H — Hrfc + He 



(5.1) 



where 



h N 



a; n 



a=l 71=1 

h N-l 

+ v\a] n) {a; n + 1\ + |a; n){a;n — 1\ 



h N 



(5.2) 



a=l n=l 



a=l n=2 



and 



He = X^tat^(|0)(a;l| + |a;l)(0|) 



(5.3) 



a=l 
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5.1. DISORDER 



It is clear from Eqs.(5.2) and (5.3) that H^c gives us the tight-binding model 
for h decoupled chains of sites (plus the origin) whereas He represents the 
coupling between the origin and the first site of each chain. 
Up to now, we set the level energy in each site (also the origin) equal to 
^0 = = and no disorder was present in this system. Here, we want to 
introduce the disorder in such a closed system, in order to apply the effec- 
tive Hamiltonian technique to the case of random variations in the diagonal 
energie^ ± 5Eq, where 5Eq is a random variable uniformly distributed in 
\—W/2, -\-W/2\ and is a disorder parameter. For a one- dimensional chain, 
there is a theoretical model, named Anderson model that leads to Anderson 
localization. Anderson localization, also known as strong localization, is the 
absence of diffusion of waves in a disordered medium. This phenomenon is 
named after the American physicist P. W. Anderson, who was the first one 
to suggest the possibility of electron localization inside a semiconductor, pro- 
vided that the degree of randomness of the impurities or defects is sufficiently 
larg^ Anderson localization is a general wave phenomenon that applies to 
the transport of electromagnetic waves, acoustic waves, quantum waves, spin 
waves, etc. 

The eigenstate of the Anderson model are exponentially localized on the 
system sites, with exponential tails given by exp(— x/La„d), where, for weak 
disorder, the localization length Land at the center of the energy band can 
be written as 

L„„rf^l05.2( — j (5.4) 

For Land ^ N {N is the number of the sites in the 1-d chain) we have the 
localized regime. The condition Land = N defines a critical value of (W/v)^^ 
for the localized regime, at any given A^. Obviously, our system is not a 



1-d chain; however, the equation (5.4) can provide a rough estimate of the 
critical value of disorder in which the system moves to the localized regime. 
We define the wavefunction of the star graph system as the eigenfunction 



-"^Also in this chapter, we set Eq = Sa = 0, but here SEq is different from zero. 
2See Rei.m. 
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5.1. DISORDER 



of the stationary Schrodinger: 

HIV'.) =^#.) (5.5) 



where H is the same defined in Eq.(5.1) and Ei are the energy levels of the 



system. Writing the eigenstates iV'j) on the site basi^ {\j)}> J = 0, . . . , hN 
we obtain 

hN 
j=0 

We define as a measure of localization the Inverse Partecipation Ratio (I.P.R.) 
of the state which is 



The value of I.P.R.* can be understood as a measure of the localization of 

the state In the limiting case of complete delocalization of state i we 

, -1 




have = (^y/hN + l) Vj that gives us I.P.R.* = 1 Vz. In the opposite 



case (complete localization) we have I.P.R.* = {hN + 1)~ In Fig.(5.1) 
we plot I.P.R.* as a function of the disorder parameter W for all the states 
within Bloch band (blue circles) and for one of the two localized states outside 
Bloch band (red dries). We considered the 5*4 (30) system with = 20 Va 
and V = 1. With these parameter values, we have (121)^^ < I.P.R. < 1. We 
notice that with increasing disorder, the delocalized states within the Bloch 
band become localized, while the localized states remain localized, despite 
the disorder. 

However, we point out that this is only a preliminary result] further study 
will indicate whether this is a general behavior of our system or not and how 
the system behavior varies according to the parameters t^, v and A^. 



^The site basis defined here is (for convenience) different from tliat adopted so far, see 



for example Eq.(5.2) and Eq.(5.3) 
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5.2. FUTURE DEVELOPMENTS 



IPR as a function of disorder W 
8^(30) system with t^= t^ = t^ = t^ = 20, v=l 




W (disorder) 

Figure 5.1: The l.P.R. as a function of the disorder parameter W for the 
5*4(30) system with ta = 20 Va and v = 1. In blue, we plot the l.P.R. of the 
states within Bloch band, [-2,2] in this case, while in green we plot one of 
the two localized outside Bloch band (the other is the same). It is evident 
the difference in the l.P.R. between the two types of states. The quantities 
plotted are dimensionless. 

5.2 Future developments 

There are many possible developments in the study of our star graph system. 
First of all, we would like to study in detail the system's behavior in the 
presence of disorder in the closed system (changing all the parameters of the 
system), but also introducing the coupling with the outside world, namely 
7 7^ 0. Secondly, we are going to study the system with asymmetric coupling 
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5.2. FUTURE DEVELOPMENTS 



with the environment and also a system with differents v from chain to 
chain. Thirdly, it would be very interesting to find possible applications of 
our model, maybe realizing an experiment that could support the theoretical 
and numerical results found. Finally, it would be of great interest the study 
of the superradiance transition in a many-body framework and also seek to 
apply the superradiance model to the photosynthesis, in order to understand 
how plants realize the coherent quantum transport at room temperature. 
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A great discovery solves a great problem, 
but there is a grain of discovery in the so- 
lution of any problem. 




George Polya - Hungarian mathematician 



Conclusions 



In this thesis, we have studied in detail the star graph system (ie chains of 
sites coupled with the origin) and we have found several unpublished results, 
both for the closed and the open system. Regarding to the open system, we 
used the effective non-Hermitian Hamiltonian approach that enabled us to 
correctly describe the coupling of the internal system to the outside world. 
Here, we divide the results obtained in two parts: closed and open system. 

Closed system 

Using the tight-binding approach, we found the eigenvalues and the eigen- 
vectors for the star graph system, also with different coupling between the 
chains and the origin. We analitically showed that there are two states out- 
side of a normal Bloch band and we also demonstrated that these two states 
are localized; we also provided the localization lengths of these states as a 
function of the parameters considered. We took in exam chains with a finite 
number of sites and also chains with a very large numbers of site. 

Open system 

Using the effective non-Hermitian Hamiltonian, we studied the eigenvalues 
and the localization lengths of the eigenstates by varying the coupling with 
the environment. We found that the localized states are very little affected 
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by the opening and we also provided an analytical estimation, for small 
coupling, of the decay rate of these states. Furthermore, we investigated 
how transport properties depend on the coupling with the outside world. 
We found an analytical estimation of the value of the coupling at which the 
superradiance transition occurs, which is also valid for chains of different 
lengths among them and with different couplings with the origin. Moreover, 
we have shown through numerical simulations that, for large number of sites 
in each chain, the maximum of the integrated transmission coincides with 
the value of the coupling at which the transition to superradiance occurs. 
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A 

Recurrence relations 



In mathematics, a recurrence relation is an equation that recursively defines a 
sequence: each term of the sequence is defined as a function of the preceding 
terms. In particular, we are interested in linear homogeneous recurrence 
relations with constant coefficients. 

An order d linear homogeneous recurrence relation with constant coefficients 
is an equation of the form: 



where the d coefficients q (for all i) are constants. 

More precisely, this is an infinite list of simultaneous linear equations, one 
for each n > d—1. A sequence which satisfies a relation of this form is called 
a linear recursive sequence or LRS. There arc d degrees of freedom for LRS, 
the initial values a^, . . . , a^^i can be taken to be any values, but then the 
linear recurrence determines the sequence uniquely. 
The same coefficients yield the characteristic polynomial 



whose d roots play a crucial role in finding and understanding the sequences 
satisfying the recurrence. 

Suppose A is a root of p{t) having multiphcity r. This is to say that {t — Xy 
divides p{t). The following two properties hold: 



On = cia„_i + C2a„_2 + . . . + caan-d 



(A.l) 



p{t) ^f- cf"-' - c^t"-^ -...-Cd 



(A.2) 
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1. Each of the r sequences A", nA", n^A", . . . , n'' ^A" satisfies the recur- 
rence relation. 

2. Any sequence satisfying the recurrence relation can be written uniquely 
as a linear combination of solutions constructed in part 1 as A varies 
over all distinct roots of 

As a result of this statement a linear homogeneous recurrence relation with 
constant coefficients can be solved in the following manner: 

1. Find the characteristic polynomial p{t). 

2. Find the roots of p{t) counting multiplicity. 

3. Write linear combination of all the roots (counting multiplic- 
ity as shown above) with unknown coefficients 6^. This is the general 
solution to the original recurrence relation. 

^ ^ (A.3) 

+ ... + (6d_,+iA: + ... + 6dn^-U:) 

where q is the multiplicity of A*. 

4. Equate each oq, Oi, . . . , from part 3 (plugging in n = 0, . . . , d into 
the general solution of the recurrence relation) with the known values 
ao, oi, . . . , from the original recurrence relation. Note, however, that 
the values a„ from the original recurrence relation used do not have 
to be contiguous, just d of them are needed (ie, for an original linear 
homogeneous recurrence relation of order 3 one could use the values 
oo, ai, 04). This process will produce a linear system of d equations 
with d unknowns. Solving these equations for the unknown coefficients 
of the general solution and plugging these values back into the general 
solution will produce the particular solution to the original recurrence 
relation that fits the original recurrence relation's initial conditions (as 
well as all subsequent values of the original recurrence relation) . 
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Interestingly, the method for solving linear differential equations is similar 
to the method above - the "educated guess" (ansai^j^ for linear differential 
equations with constant coefficients is e^^ where A is a complex number that 
is determined by substituting the guess into the differential equation. This 
is not a coincidence. If you consider the Taylor series of the solution to a 
linear differential equation: 



you see that the coefficients of the series are given by the n derivative of 
f{x) evaluated at the point a. The differential equation provides a linear 
difference equation relating these coefficients. This equivalence can be used 
to quickly solve for the recurrence relationship for the coefficients in the 
power series solution of a linear differential equation. 
Let us give a very famous example of homogenous recurrence relation. 

Example: Fibonacci sequence 

The Fibonacci sequence is named after Leonardo of Pisa, who was known as 
Fibonacci (a contraction of filius Bonacci, "son of Bonaccio"). Fibonacci's 
1202 book Liber Abaci introduced the sequence to Western European math- 
ematics, although the sequence was independently described in Indian math- 
ematics and it is disputed which came first. 

Fibonacci numbers are used in the analysis of financial markets, in strategies 
such as Fibonacci retracement, and are used in computer algorithms such as 
the Fibonacci search technique and the Fibonacci heap data structure. The 
simple recursion of Fibonacci numbers has also inspired a family of recursive 
graphs called Fibonacci cubes for interconnecting parallel and distributed 
systems. They also appear in biological settings, such as branching in trees, 
arrangement of leaves on a stem, the fruit spouts of a pineapple, the flowering 

^ In physics and mathematics, an ansatz is an educated guess that is verified later 
by its results. An ansatz is the establishment of the starting equation(s) describing a 
mathematical or physical problem. It can take into consideration boundary conditions. 




(A.4) 



n=0 
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of artichoke, an uncurling fern and the arrangement of a pine cone. 

In mathematics, the Fibonacci numbers are the numbers in the following 

integer sequence: 



0, 1, 1, 2, 3, 5, 8, 13, 21, 34, 55, 89, 144 .. . 



(A.5) 



By definition, the first two Fibonacci numbers are and 1, and each subse- 
quent number is the sum of the previous two. In mathematical terms, the 
sequence F„ of Fibonacci numbers is defined by the recurrence relation 



F — F 1 + F 



with seed values 



Fn = 



The characteristic equation is 



with roots 



a — 



1 + V5 



1 - V5 



We have 



Fn = Aa"" + 
From the initial conditions Fq = and Fi = 1, 

Fo^A + B^O 
Fi^Aa + Bp^l 

Solve for A and B, we have A = 1/^5 and B = -1/^5. 
Hence 



F — 



1 + V5 



1 - VE 



(A.6) 



(A.7) 



(A.8) 



(A.9) 



(A.IO) 



(A.ll) 



(A.12) 



2 / \ 2 

In this way, we obtained a closed-form solution, known as Binet 's formula. 
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B 

Determinant of a tridiagonal matrix 



A tridiagonal matrix is a matrix that has nonzero elements only in the main 
diagonal, the first diagonal below this, and the first diagonal above the main 
diagonal. This is the most general tridiagonal matrix: 



M 





ai hi 





\ 




Ci ■•• 









■•• 


' ■ • hn-1 




\ 





Cn— 1 Oifi 


/ 



(B.l) 



For tridiagonal matrices, determinants can be evaluated via multiplication 
of 2 X 2 matrices as follow^ 

f ai bi \ 



det Ta 



M 




(B.2) 




11 



where the subscript 11 indicates that we have to consider the first element 
in the first column. 



iSee Ref.IMl- 
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